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I.  INTRODUCTION 


In  this  study,  the  lift  and  drag  coefficients  are  calculated  as 
a  function  of  the  ratio  of  the  cylinder's  tangential  velocity  to  the 
velocity  of  the  free  streaai  flow.  Here  the  boundary- layer  is  laminar 
to  separation  with  the  free  shear  layers  near  the  separation  points 
also  being  laminar.  According  to  Swanson^,  the  nondimensionalized  wall 
velocity  (peripheral  velocity  of  cylinder  divided  by  freestream  velocity) 
must  always  be  less  than  0.5  or  the  boundary  layer  will  not  be  laminar 
to  separation.  From  Wu's2  review  article  on  wakes  and  cavities,  laminar 
flow  for  the  shear  layers  near  the  separation  points  requires  flows  with 
Reynolds  numbers  less  than  5*10  .  Yet  the  Reynolds  number  must  be  high 
enough  for  the  boundary- layer  approximation  to  be  valid;  that  is,  the 
boundary- layer  is  very  thin  and  the  pressure  immediately  outside  the 
boundary- layer  is  impressed  through  the  boundary  layer  thickness. 

Rotating  a  body  of  revolution  so  that  its  axis  of  rotation  is  at 
an  angle  with  its  direction  of  relative  fluid  velocity  not  only  affects 
the  drag  but  also  introduces  a  new  force.  This  force,  called  the  Magnus 
force,  is  directed  perpendicular  to  the  plane  in  which  the  rotational 
ixis  and  direction  of  translation  lie.  G.  Magnus^  correctly  attributed 
this  force  to  the  pressure  field  produced  by  the  inviscid-velocity  dis¬ 
tribution  about  the  rotating  body  in  the  airstream.  He  did  this,  and 
set  the  tone  for  subsequent  experiments,  by  considering  the  apparently 
simpler  two-dimensional  problem  of  the  rotating  cylinder  in  crossflow. 
Using  this  model,  he  established  that  the  force  was  directed  toward 
the  side  of  the  cylinder  moving  with  the  direction  of  flow. 

4 

Lord  Rayleigh  was  the  first  to  construct  a  mathematical  model  of 
the  flow  field  by  assuming  an  ideal  fluid;  a  potential  vortex  com¬ 
bined  with  a  doublet  in  a  uniform  flow  induces  a  circulation  and  pro- 


1.  W.  M.  Swanson,  "An  Experimental  Investigation  of  the  Two- 
Dimensional  Magnus  Effect , "  Final  Report ,  Office  of  Ordnance 
Research  Contract  No.  DA-33-019-ORD-1434 ,  Dept,  of  Mech  Eng., 

Case  Institute  of  Technology ,  33  De  -'ember  1956. 

2.  T.  Y,  Wu,  "Cavity  and  Wake  Flows,"  in  Annual  Review  of  Fluid 
Mechanics .  Vol.  IV,  Van  Dyke,  Vincenti,  and  Weliausen  (eds.). 

Annual  Reviews ,  Inc.,  Palo  Alto,  Calif.,  1972,  pp.  243-284. 

3.  G.  Magnus,  "On  the  Deflection  of  a  Projectile,"  English 
translation  in  Taylor's  Scientific  Memoirs,  1853. 

4.  Lord  J.  W.  S.  Rayleigh ,  "On  the  Irregular  Flight  of  a  Tennis  Ball," 
Scientific  Papers  1,  1857 ,  pp.  344-346. 
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duces  a  side  force  according  to  Zhukovskii's  theorem.  However,  no 
mathematical  method  could  yet  be  used  to  couple  the  rotation  of  the 
cylinder  with  the  fluid  to  produce  the  circulation.  A  missing  element 
was  found  when  Prandtl5  introduced  the  concept  of  the  boundary- layer 
in  1904,  thus  providing  a  better  understanding  of  the  rotating 
cylinder  problem,  Prandtl  and  Tietjens6  also  showed  that  separation 
was  delayed  on  the  upper  part  of  the  cylinder  wall  moving  with  the 
flow  while  separation  was  hastened  on  the  underside  of  the  cylinder. 
However,  no  successful  theoretical  solutions  were  found  for  quanti¬ 
fying  the  lift  and  drag  forces. 

.7 

Reid  measured  the  lift  and  drag  coefficients  of  the  rotating 

cylinder  for  rotational  velocities  up  to  u  =  3.4,  where  u  is  the 

w  w 

tangential  velocity  of  the  surface  of  the  cylinder  divided  by  the 
free  stream  velocity.  A  maximum  lift  to  drag  ratio  of  7.8  was  observed, 
and  it  was  noted  that  the  drag  decreased  somewhat  with  increasing 
velocity  of  rotation.  A.  Thora®»®»^®»H»  12,13  and  hjs  associates  math 
comprehensive  studies  involving  the  variations  of  the  Magnus  force 
with  Reynolds  number,  surface  condition,  end  conditions  and  other 
factors.  Pressure  field  data  were  also  obtained  with  the  lift  and 


S.  r..  Prandtl,  " Uber  Flussigkeitsbewegung  bei  sekr  Weiner  Rcilur.g," 
Verb.  III.  int.  Math  Kongr .,  Heidelberg,  1904,  pp.  484-421. 
Bnglieh  translation  available  as  NACA  TM-452. 

f Prandtl  and  0.  G.  Tietjens,  Applied  Hydro-  and  Aeromechanics, 
New  York,  Dover  Publications,  Inc.,  1957 ,  p.‘  82~  ' - 

'•  E’  “•  I'bid,  "Tests  of  Rotating  Cylinders,"  NACA  TN  £09,  1224. 

8.  /!.  Thom,  "The  Aerodynamics  of  a  Rotating  Cy Zuider,"  Ph.V. 
Dissertation ,  University  of  Glasgow,  19*26.' 

Thom,  "Experiments  on  the  Air  Forces  on  Rotating  Cylinders 
ARC  R4I1  101 P,  1925. 

10.  A.  Thom,  "The  Pressures  Round  a  Cylinder  Rotatim  in  an  Air 
Current,"  ARC  RAM  1082,  1926. 

11.  A.  Thorn,  "Experiments  on  the  Flow  Past  a  Rotating  Cylinder 
ARC  RAM  1410,  1951. 

A.  1  horn  and  S.  R.  Scngupta,  "Air  Torque  on  a  Cylinder  Rotating  in 
an  Air  Stream,  "  ARC  RAM  1520,  1932. 

13.  A.  Thom,  " Effects  of  Discs  on  the  Air  Forces  on  a  Rotatim 
Cylinder,  "  ARC  RAM  1623,  1934. 
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drag  coefficients  being  deduced  from  these.  More  recently,  Swanson 
has  investigated  the  rotating  cylinder  and  has  presented  experimental 
results  for  a  wideT  range  of  rotation  rates  and  Reynolds  numbers  than 
were  previously  available.  In  his  careful  study,  three-dimensional 
effects  due  to  finite  length  of  the  cylinder  were  minimized  and 
pressure  probe  results  for  the  boundary- layer  were  obtained  in  addition 
to  lift  and  drag. 

Figure  1,  obtained  from  Swanson's  work,  illustrates  some  features 
of  the  flow  in  the  boundary  layer  for  a  wall  velocity  to  free-stream 
velocity  ratio  (u^)  of  one.  Caution  must  be  used  in  applying  some 

features  of  this  boundary-layer  flow  to  the  problem  under  investi¬ 
gation  since  Swanson  concluded  that  the  boundary  layer  on  the  under¬ 
side  of  the  cylinder  for  this  rotation  rate  was  very  probably  turbulent 
before  separation.  The  stagnation  point  location  corresponding  to  this 
velocity  ratio  is  nearly  the  same  as  the  location  for  the  nonrotating 
cylinder.  On  the  top  of  the  cylinder,  at  an  angle  of  about  120°  from 
the  stagnation  point,  the  boundary  layer  thickens  rapidly  and  separation 
appears  to  occur  with  a  velocity  profile  satisfying  the  Moore-1  Rott  - 
Sears  (unpublished)  criterion  of  u  =  =  0  at  a  coincident  point.  A 

boundary  layer  then  appears  and  grows  on  the  cylinder  in  the  wake  start¬ 
ing  at  the  upper  separation  point.  The  profile  corresponding  to  the 
apparent  lower  separation  point  does  not  obey  the  Moore-Rott -Sears 
criterion  for  separation;  from  the  form  of  the  profiles,  it  is  not 
obvious  what  the  correct  criterion  would  he.  But,  as  noted  before,  the 
boundary  layer  is  probably  turbulent  in  this  region. 

Griffiths  and  Ma16  have  investigate  ’•.he  Magnus  force  phenomena 
particularly  with  regard  to  the  negative  Magnus  force  and  its  rela¬ 
tionship  to  the  Reynolds  number.  The  negative  force  is  thought  to 
occur  because  the  side  of  the  cylinder  going  against  the  flow  will 
have  a  higher  local  relative  Reynolds  number  for  corresponding  points 
than  the  other  side.  Here  the  relative  Reynolds  number  is  defined  as 


14.  F.  K.  Moore ,  "Chi  the  Separation  of  the  Unsteady  Laminar  Boundary 
Layer, "  Proceedings  of  Symposium  on  Boundary  Layer  Research  held 
on  August  1957,  H.  Gdrtler  (ed.),  Spring er-Verlag  Press,  Berlin, 
1958. 

15.  N.  Rott  with  F.  K.  Moore,  (ed.)  Theory  of  Laminar  Flows, 
Princeton  University  Press,  Princeton,  R.J.,  1964,  p.  452. 

16.  R.  T.  Griffiths  and  C.  Y.  Ma,  "Differential  Boundary-Layer 
Separation  Effects  in  the  Flow  over  a  Rotating  Cylinder, " 

Aero.  J.  of  the  Roy.  Aero.  Soc .,  Vol.  73,  1969,  pp.  521-526. 
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Figure  1.  Boundary-Layer  Profiles  on  a  Rotating  Cylinder 
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uq  x  (1  ?  uw/uo)/v,  where  the  minus  sign  corresponds  to  the  upper  side 

and  the  plus  sign  refers  to  the  lower  side.  If  the  Reynolds  number  for 
the  flow  is  in  a  certain  range,  transition  to  turbulence  will  occur  on 
the  upstream  moving  wall  thus  delaying  separation  and  causing  the  cir¬ 
culation  strength  to  be  of  opposite  sign.  The  present  theoretical 
investigation,  however,  is  not  concerned  with  the  turbulent  boundary 
layer  and  this  interesting  aspect  of  Magnus  forces  will  not  be  discussed 
further.  Nevertheless,  Griffiths  and  Ma's  results  would  seem  suspect 
in  a  quantitative  sense  because  their  lift  and  drag  forces  do  not 
agree  for  common  regions  of  the  Reynolds  number  with  those  of  Thom 
and  Swanson's.  In  fact,  their  drag  coefficient  for  the  non-rotating 
cylinder  is  some  30%  lower  than  found  by  other  investigators.  Conse¬ 
quently,  their  results  will  not  be  used  for  comparison  purposes  in 
this  study. 

The  most  successful  attempt  to  theoretically  treat  the  Magnus 
force  (for  any  rotation  rate)  is  M.  B.  Glauert's*7  treatment  of  a 
cylinder  whose  peripheral  velocity  is  large  enough  so  that  the  stream¬ 
lines  around  the  immediate  vicinity  of  the  cylinder  are  closed.  He 
then  assumes  that  the  velocity  outside  the  boundary  layer  is 

fl  =  fl  +  2u  sin  x, 
e  c  o  ’ 


where  u  is  the  free  stream  speed,  flc  is  the  circulatory  component 

of  the  fluid  velocity  and  x  is  the  angular  displacement  from  the  line 
through  the  center  of  the  cylinder  aligned  with  the  direction  of  the 
unperturbed  flow.  He  further  assumes  a  parameter  perturbation  ex¬ 
pansion  together  with  an  expansion  in  terms  of  exp(ix)  and  substitutes 
the  expansion  into  the  boundary  layer  equations.  He  solves  the  equations 
and  obtains 


2tt  a  u 


w 


.1-3  (u0/uu)2  -  3.24  (u/u/ 


where  r  is  the  circulation  strength  and  a  is  the  radius  of  the 
cylinder.  This  indicates  that  the  lift  force  =  p  uQr  is  asymptotical¬ 
ly  a  linear  function  of  the  peripheral  velocity  of  the  cylinder  and  has 
no  upper  limit.  This  result  contradicts  Prandtl's*®  assertion  that  the 


17.  U.  B.  Glauert,  "The  Flow  Past  a  Rapidly  Rotating  Circular  Cylinder," 
Royal  Soo.  of  London  Proceedings,  Vol.  A-242,  1957,  pp.  103-115. 

18.  L.  Prandtl,  Die  Naturwissenschaften,  Vol.  IS,  1925,  pp.  93-108. 
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upper  limit  for  the  lift  coefficient  is  4.  Glauert  is  supported  by 
Swanson’s  experimental  results;  Prandtl's  upper  limit  is  exceeded  and 
no  other  upper  limit  was  encountered. 

In  contrast,  some  other  attempts  to  predict  the  lift  and  drag 
coefficients  need  data  from  experiment  to  obtain  meaningful  results. 

W.  G,  Bickley1^  placed  a  single  vortex  downstream  of  a  cylinder  with 
circulation  superimposed  upon  a  flow  around  a  cylinder.  The  location 
of  this  vortex  was  selected  to  give  the  experimentally  obtained  rela¬ 
tionship  between  CL  and  Cp  for  large  values  of  C^.  T.  Gustafson20 

expanded  Bickley's  analysis,  but  distributed  the  vorticity  on  the 
separation  streamlines  extending  to  infinity  instead  of  using  a  single 
concentrated  vortex.  The  vortices  traveling  downstream  result  in  net 
drag  and  lift  forces  on  the  cylinder.  Gustafson  was  also  concerned 
with  modeling  the  pressure  distribution  upstream  of  separation  but 
Bickley  was  not  concerned  with  such  details  of  the  flow  field. 

The  inviscid-flow  model  used  in  the  present  investigation  is  a 
specialized  case  of  Piercy,  Preston  and  Whitehead's^  model.  Their 
model  was  used  to  consider  the  flow  around  other  shapes  such  as 
ellipses  and  flat  plates  set  at  various  angles  to  the  direction  of 
free  stream  flow.  Two  complex  planes  are  considered:  the  physical 
plane  and  the  untransformed  plane.  The  inviscid  flow  in  the  untrans¬ 
formed  plane,  shown  in  Figure  2,  results  from  the  following  super¬ 
positions  onto  a  uniform  flow:  a  doublet,  two  sources  located  on  the 
downstream  part  of  the  cylinder  and  also  rear  the  rear  of  the  cylinder, 
their  image  sinks  at  the  center  of  the  cylinder,  and  a  point  vortex 
located  in  the  center.  With  this  configuration,  the  separation  stream¬ 
lines  emanate  from  the  cylinder  at  right  angles  to  the  surface.  This 
complex  potential  field  is  conformally  mapped  onto  the  physical  plane 
by  a  Zhukovskii  transformation  that  doubles  the  angle  of  the  stream¬ 
lines  emanating  from  the  cylinder  at  Sj  and  S2.  The  cylinder  in  the 


19.  W.  G.  Bickley,  " The  Influence  of  Vortices  Upon  the  Resistance 
Experienced  by  Solids  Moving  Through  a  Liquid , "  Proc,  Roy. 

Soa.  of  London,  Vol.  A-119 ,  1928 ,  pp.  146-156. 

20.  T.  Gustafson ,  "On  the  Magnus  Effect  According  to  the  Asymptotic 
Hydrodynamic  Theory , "  Edkan  Ohlssons  Buchdruckerei ,  Lund ,  Sweden, 
1933.  NACA  translation  N-25921 ,  1954. 

21.  N.  A.  Piercy ,  J.  H.  Preston,  and  L.  G.  Whitehead ,  "The  Approximate 
Prediction  of  Skin-Friction  and  Lift,"  Phil,  Mag.  Vol.  26,  1938, 
pp.  791-815. 
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i  the  Untransfomed  Planes  (Source-Wake 

) 


5 


physical  plane  together  with  the  corresponding  streamlines  are  shown 
in  Figure  3„  For  purposes  of  calculating  the  drag  and  the  lift  forces, 
the  flow  inside  the  separation  streamlines  is  ignored  and  the  base 
pressure  is  assumed  constant  at  the  separation  value.  The  pressure 
gradients  near  the  separation  points  are  required  to  be  finite  as  is 
found  experimentally.  In  this  inviscid  flow  model,  the  pressure  values 
on  the  separation  streamlines  increase  asymptotically  towards  the  free- 
stream  value  as  actually  occurs.  The  separation  streamlines  asymp¬ 
totically  approach  a  finite  width  apart  with  distance  from  the  cylinder. 
It  is  assumed  that  this  inviscid  flow  model  simulates  the  shear-layers 
near  the  separation  points  and  hence  produces  a  more  realistic  velocity 
distribution  on  the  cylinder  ahead  of  separation. 

The  general  approach  defined  above  was  used  for  the  problem  of 
the  nonrotating  cylinder  in  crossflow  by  Bluston  and  Paulson22.  They 
employed  the  source-wake  of  Parkinson  and  Jandali23>  a  specialized 
model  of  the  Piercy,  Preston  and  Whitehead  model.  Their  inviscid-f low 
model  in  the  untransformed  plane  had  t’.-?o  equal  sources  located  sym¬ 
metrically  about  the  flow  axis  together  with  their  image  sinks  at  the 
center  of  the  circle  to  produce  symmetric  separation  streamlines.  The 
resultant  complex  potential  field  was  transformed  to  the  physical  plane 
by  a  Zhukovskii  transformation.  Parkinson  and  Jandali  needed  the 
velocity  at  separation  and  the  location  of  the  separation  point  to 
specify  their  inviscid  flow  model.  However,  Bluston  and  Paulson  ob¬ 
tained  their  inviscid-flow  model  by  specifying  the  separation  point 
and  requiring  that  the  pressure  gradient  be  finite  at  separation. 

The  inviscid-flow  model  used  by  Bluston  and  Paulson  was  made  to 
agree  with  the  results  of  boundary- layer  theory.  Using  the  resulting 
velocity  distribution  for  a  trial  value  of  the  inviscid-flow  separation 
point,  they  calculated  the  boundary  layer  over  the  surface  of  the 
cylinder  until  boundary  layer  separation  was  reached.  A  second  in¬ 
viscid  velocity  distribution  was  generated  with  the  inviscid-flow 
separation  point  determined  from  the  first  boundary -layer  calculation. 
This  procedure  was  repeated  until  the  location  of  the  separation  points 
calculated  by  boundary  layer  theory  agreed  approximately  with  the 
separation  points  of  the  inviscid  flow  model.  Separation  was  achieved 
at  83°  compared  with  about  82°  for  experiment;  furthermore,  the  pressure 
distribution  over  the  cylinder  ahead  of  separation  was  in  good  agree¬ 
ment  with  experimental  measurements.  These  results  suggested  the 


22.  H.  S.  Bluston  and  R.  W.  Paulson ,  nA  Theoretical  Solution  for 
Laminar  Flow  Past  a  Bluff  Body  with  a  Separated  Wake ,  n  Journal 
de  Mecanigue ,  Vol.  2 ,  No „  1 ,  1972 ,  pp.  161-180.  ~~~  “ 

28.  G.  V.  Parkinson  and  T.  Jandali,  "A  Wake  Source  Model  for  Bluff 
Body  Potential  Flow,"  J.  Fluid  Mech.  Vol.  40.  Part  3.  1970. 

pp.  97 7-594.  '  "  ’  ~~ 
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Figure  3.  Flow  Past  Rotating  Cylinder  in  the  Physical  Planes 
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possibility  of  extending  the  technique  to  the  more  difficult  problem  of  the 
rotating  cylinder  in  crossflow. 

Although  the  approach  of  the  present  problem  is  similar  to  Bluston 
and  Paulson's  approach,  the  problem  addressed  in  the  present  report  is 
more  difficult  then  their  problem  in  at  least  two  ways.  The  inviscid- 
f low  model  used  is  complicated  by  the  loss  of  symmetry  (unequal  source 
strengths  located  asymmetrically  with  respect  to  the  free  stream  flow) 
and  the  addition  of  a  point  vortex.  A  computer  had  to  be  used  to  find 
five  unknown  parameters  in  the  untransformed  plane  whereas,  the  param¬ 
eters  for  Bluston  and  Paulson's  problem  were  found  analytically.  The 
other  difficulty  involved  calculating  the  boundary  layer  on  the  rotating 
cylinder.  For  the  part  of  the  cylinder  where  the  wall  is  moving  against 
the  flow,  there  is  a  region  in  the  boundary  layer  near  the  wall  which 
has  reverse  flow  and  calculations  of  boundary  layers  with  regions  of 
reverse  flow  have  seldom  been  attempted.  However,  an  integral  technique 
has  been  used  to  numerically  compute  the  boundary  layer  for  such  a 
situation24  in  which  the  integral -momentum  equation  and  the  integral- 
energy  equation  were  s  multaneously  integrated.  The  moving-wall  simi¬ 
larity  solutions25  were  used  to  approximate  the  nonsimilar  velocity  pro¬ 
files  encountered  on  the  rotating  cylinder. 


IX.  THE  INVISCID-FLOW  MODEL 
A.  General  Description 

The  two-dimensional,  incompressible,  irrotational  point  vortex 
source-wake  model  for  a  spinning  cylinder  is  represented  in  the  physical 
plane  in  Figure  3  where  ug  is  the  dimensionless  inviscid  velocity  tangen¬ 
tial  to  the  surface,  uw  is  the  peripheral  dimensionless  velocity  of  the 
cylinder  rotating  in  the  clockwise  sense,  uq  is  the  free  stream  velocity 

that  is  used  to  nondimensionalize  Q  and  fl  ,  and  x  is  a  nondimensionalized 

e  w 

coordinate  which  defines  an  angular  position  on  the  cylinder.  The 
cylinder  and  general  flow  field  are  described  in  two  complex  planes: 
the  z  complex  plane  having  its  real  axis  aligned  with  the  direction  of 
unperturbed  flow  and  the  z  plane  having  its  imaginary  axis  going  through 


24.  K.  S.  Pansier  and  J.  E.  Danberg,  "An  Integral  Analysis  of  Boundary. 
Layers  on  Moving  Vails,"  U.  S.  Amy  Ballistics  Besearch  Laboratories 
Report  1792,  June  1975.  AD  M01220S 

55.  J.  E.  Danberg  and  K.  S .  Fancier,  "Similarity  Solutions  of  the 
Boundary  Layer  Equations  for  Flow  over  a  Moving  Wall,"  U.S. 

Amy  Ballistic  Research  Laboratories  Report  1714,  April  1974, 

AD  781496. 
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the  points  of  breakaway,  Sj  and  S2.  The  center  of  the  cylinder  is  a 

point  on  the  real  axis  of  the  z  plane.  The  z  plane  is  used  in  an 
intermediate  step  of  the  transformation  from  the  unfcransformed  plane 
C  to  the  description  of  the  flow  in  the  z  plane.  At  the  points  and 

S^,  the  free  streamlines  leave  the  cylinder  smoothly  and  are  a  finite 

distance  apart  at  infinity.  The  stagnation  point  will  not  generally 
be  at  x  a  0,  and  the  flow  field  will  normally  be  asymmetrical  with 
respect  to  the  axis. 

To  put  this  source-wake  model  in  better  perspective,  wake  flow 
for  nonrotating  cylinders  with  Reynolds  numbers  between  1500  and  10 
will  be  discussed  briefly.  This  is  a  subcritical  region  within  which 
the  boundary  layer  on  the  cylinder  is  nowhere  turbulent.  According  to 
Roshko  and  Fiszdon26>  shear  layers  separate  from  the  body  and  envelop 
a  region  of  recirculating  flow  called  the  near  wake.  These  free  or 
separating  shear  layers  are  thin  and  well  defined  near  the  cylinder; 
the  flow  outside  the  free  shear  layers  is  irrotational.  The  near  wake 
is  unstable  and  oscillates  periodically,  particularly  close  to  the  end 
of  the  near  wake  region;  the  shear  layers  roll  up  into  large  vortices 
near  the  closure  of  the  near-wake  and  then  break  away  creating  a 
Kantian  vortex  street,  Behind  this  near  wake,  the  flow,  being  vortical 
and  turbulent,  is  called  the  turbulent  far-wake.  As  the  Reynolds 
number  is  increased,  the  tr-  .sition  to  turbulence  advances  along  the 
free  shear  layers  towards  the  points  of  separation.  The  time-averaged 
pressure  is  found  to  have  an  almost  constant  value  over  the  part  of 
the  cylinder  immersed  in  the  wake. 

The  source-wake  bound-vortex  model  used  in  this  investigation 
does  take  into  account  some  phenomena  observed.  The  free  streamlines 
of  the  model  take  the  place  of  the  time-averaged  free-shear  layers 
and  form  boundaries  of  the  irrotational  flow  external  to  the  wake 
region.  Also,  the  fluid  velocities  on  the  free  stream-lines  decrease 
with  distance  from  the  cylinder  towards  the  free  stream  velocity  as 
observed  in  experiments.  It  is  found,  however,  that  streamlines 
continually  enter  the  actual  wake  which  widens  downstream  because  of 
the  diffusion  of  vorticity.  Also  the  near  wake  is  continuously 
shedding  Karnian  vortex  streets  at  the  rear  of  the  bubble  and  thus 
the  free  shear  layers  are  unsteady  over  part  of  the  near-wake  region. 

A  successful  treatment  to  take  into  account  the  diffusion  of  vorticity 
and  unsteadiness  in  the  wake  has  not  been  developed.  A  flow  model 
taking  into  account  these  features  of  the  wake  would  require  an  ex¬ 
tensive  investigation  of  the  wake  region.  Since  the  intention  is 


26.  A.  Boshho  and  W.  Fiszdon ,  Problems  of  Hydrodynamics  and  Continum 
Meohanio8t  Soo.  Ind.  Appl.  Math ,  Philadelphia ,  1967,  pp.  606-616. 
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not  to  study  the  wake  itself  but  rather  to  account  for  the  effect  of 
the  wake  on  the  upstream  (attached  flow)  pressure  distribution.,  it  is 
expected  that  a  reasonable  model  for  the  time-averaged  near-wake  free 
shear  layers  would  result  in  an  improvement  over  their  complete  neglect. 

To  construct  the  flow  field  just  described,  flow  in  an  untransformed 
plane  c  is  considered  as  in  Figure  2.  The  fundamental  flow  past  the 
circle  is  the  uniform  flow  plus  a  doublet.  Additionally,  sources  of 
strength  2^  and  2Q,  arc  located  at  angles  y  and  6  respectively  from 

the  direction  of  unperturbed  flow.  Image  sinks  -Q^  and  -C^  are  placed 

at  the  origin  to  satisfy  the  boundary  conditions  on  the  circle.  A  point 
vortex  of  strength,  r,  is  also  placed  at  the  origin  and  does  not  affect 
the  boundary  conditions  on  the  circle  but  does  affect  the  positions  of 
Sj  and  S,  and  the  complex  potential  field.  Here  R  is  the  radius  of  the 

cylinder.  Now  half  the  distance  between  the  two  sources  will  be  assign¬ 
ed  the  value  unity.  Therefore,  R  is  equal  to  esc  a. 

The  complex  potential  of  the  resulting  flow  referred  to  the  c,  plane 
is : 

X(0  =  Vq  (r,  +  R2/ c)  +{ 2Qj  (In  (r.  -  Rely)  -  h  In  (r.)]  /  2n}+ 

|  2Q?  [In  (c  -  Re1<S)  -  h  In  (01  /2n}+  ( h  i  r/n)  In  ?  (1) 


For  the  c,  plane,  only  the  fluid  dynamical  properties  on  the  circle 
C  will  be  of  primary  interest.  The  velocity  on  the  circle  C  can  be 
obtained  as  shown  below.  From  Figure  2,  the  directed  vectors  from  the 
sources  2Q  and  2Q  to  the  point  x,  can  be  seen  to  be  respectively, 


„  i6  it  o'\ 

5  -  Re  =  r2e  ,  [2) 

where  s  and  t  are  angles  between  the  lines  drawn  from  the  sources  to  a 
point  r,  on  the  circle  and  the  lines  drawn  through  the  sources  parallel 
to  the  real  axis.  Here,  Tj  and  r,  are  the  lengths  from  the  point  r, 

to  the  sources  2Q,  and  2Q_  respectively.  Thus,  using  equation  (2)  to 

rewrite  the  expression  for  the  complex  potential  m  equation  (1),  the 
velocity  potential  can  then  be  expressed  as 


t> 


2  V  esc  a  cos 
o 


+  + 


—  In  (r.  /esc  a) 

2  it  A 


2ir 


(3) 
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Now  i^2  and  r*  can  be  seen  to  be  the  following  in  terms  of  esc  a,  y, 
6,  and  y: 


=  2  (csc2a)  [  1  -  cos  (<t>-y)]  , 

r,2  =  2  (csc2a)  [  1  -  ('os  (<M)] 

The  tangential  velocity  at  the  surface  of  the  cylinder  is  just 

d$ 

A  s  _  1$  sin  0 

Thus,  substituting  equation  (4)  into  equation  (3),  differentiating 
according  to  equation  (S') ,  and  using  trigonometric  identities,  the 
following  equation  can  be  obtained: 


a  =  2  V  sin 
M  o 


2ttcsc  a 


_  [Qx  cot  (4^0  +  $2  cot  (^-6-)]  +  r /2ti  esc  5.(6) 


By  using  the  following  definitions 

q  5  V\  . 

«1,2  $  5j,2  '  2'V0  “C  S  ’  (?) 

u  =  17  2vV  CSC  a  , 

C  0 

the  following  equation  can  be  obtained  from  equation  (6). 

q  =  2  sin  4-  [Qj  cot  O^— )  +  Q2  cot  +  uc  *  ^ 

Equation  (8)  is  the  expression  for  the  dimensionless  tangential  velocity 
at  a  point  on  the  surface  of  the  cylinder  in  the  untransformed  plane 

The  resulting  complex  potential  and  the  complete  circle  E  in  the 
?  plane  can  be  mapped  to  the  physical  plane  z  by  first  transforming  to 
the  l  plane.  The  complex  plane  c  has  its  real  axis  V?e 

direction  of  the  freestream  flow.  The  complex  plane  5  has  lts  reai 
axis  at  an  angle  -  e  from  the  real  axis  of' the  5  plan*  and  is  located 
such  that  S1  and  S2  are  symmetric  about  the  real  axis.  An  angle  <f 

describing  the  position  on  the  circle  in  the  c  plane  can  be  described 
in  the  l  plane  by  the  relationship 

a  =  ♦  e 
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The  r,  and  \  planes  will  be  termed  the  untrans formed  planes.  Points  on 
the  l  plane  can  be  mapped  conformally  into  points  on  the  z  plane  by  an 
analytic  function: 


z  =  M(£)  .  (9) 

The  particular  mapping  function  is  a  Zhukovskii  transformation: 


M(4) 


cot  a  - 


1 

l-  cot  a. 


(10) 


Here  cot  a  is  just  the  distance  from  the  origin  along  {  to  the  straight 
line  drawn  from  to  in  Figure  2.  The  points  S^,  S2,  are  critical 

points  and  M*  (?)  thederivative  of  M  is  zero  at_the  points.  Angles  of 
intersection  in  the  \  plane  are  doubled  in  the  z  plane  by  the  mapping 
function  at  and  S^.  Thus,  the  angle  of  intersections  of  the  separa¬ 
tion  streamlines  with  the  circle  E  in  the  ;  plane  increase  from  90  to 
a  tangential  intersection  with  the  arc  C.  Also,  the  circle  E  is  mapped 
onto  the  slit  S^AS^BS^.  The  points  and  lie  on  the  axis  of  the 

z  plane. 


The  relationship  between  the  two  complex  velocities  in  the  two 
planes  z  and  ?  is  given  by 


W(z)  =  w(c)  /  M'(0 


(11) 


where 


M’(0  =  1  +  1/U  -  cot  a)2  .  (12) 

From  equation  (11)  and  equation  (12),  it  is  apparent  that  at  large 
distances  from  the  flow,  the  complex  velocities  in  the  two  planes 
asymptotically  approach  each  other;  and  this  implies  that  uq  =  Vq  and 

e  =  e . 


The  transformation  from  the  z  plane  to  the  z  plane  is  accomplished 
by  a  translation  and  then  a  rotation  through  an  angle  e.  The  z  and  z 
planes  will  be  teimed  the  physical  planes.  As  mentioned  before,  only 
the  dynamical  quantities  on  the  cylinder  are  of  primary  interest;  hence, 
only  the  circle  E  of  Figure  2  will  be  mapped  into  the  arc  C  of  Figure  3. 

In  the  sections  that  follow,  it  will  be  necessary  to  know  the 
tangential  velocity  u?  on  the  arc  C  of  the  cylinder  in  the  z  plane  in 

terms  of  the  position  4>  on  the  circle  E.  The  tangential  velocity  on 
the  cylinder  in  the  physical  plane  z,  given  the  corresponding  velocity 
in  the  plane  ?,  can  be  obtained  if  equation  (11)  is  used  specialized  to 
the  circle: 
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(13) 


un(x)  ■  q($)/N($)  , 

6 


when  N  (4>)  is  defined  as 


N(*)=|M’(t)|  E  (14) 

The  subscript  E  means  that  M'(t)  is  evaluated  on  the  cylinder  E.  The 
value  N(<j>)  is  obtained  by  multiplying  M'(t)  on  the  circle  E  by  its  com¬ 
plex  conjugate  and  using  trigonometric  identities  to  obtain  squared 
quantities  in  both  numerator  and  denominator.  Then  taking  the  square 
root  and  putting  c  in  terms  of  *  and  e,  the  expression  for  becomes 


Ntt)  .  2  I  c°s  V  cos  - 5^  (IS) 

1-2  cos  a  cos  ($+e)  •*  cos  a 

The  inviscid  flow  velocity  can  then  be  found  by  using  equation  (13). 
The  function  for  N($),  equation  (15),  is  known  (since  a  and  e  can  be 
obtained  from  equation  (33) . 

B.  Curvature  and  Pressure  Conditions  at  the  Separation  Points 

Parkinson  and  Jandali23,  with  their  symmetric  flow  model,  related 
the  location  and  strength  of  the  source  in  the  untransformed  plane  to 
two  parameters  in  the  physical  plane:  the  position  of  separation  and 
the  'ocal  external  velocity  at  separation.  By  using  the  finite 
curvature  condition  at  the  separation  point  originally  suggested  and 
used  by  Woods27,  Bluston  and  Paulson22  could  characterize  the  inviscid 
flow  in  the  physical  plane  with  one  parameter:  the  position  of  sep¬ 
aration.  The  number  of  parameters  of  the  physical  plane  needed  to 
describe  the  present  flow  model  will  also  be  reduced  by  imposing  sim¬ 
ilar  conditions.  The  streamline  curvature  condition  imposed  by  Bluston 
and  Paulson  is  applied  to  both  separation  points  in  the  present  case 
and  thus  reduces  the  number  of  descriptive  parameters  by  two.  After 
applying  a  second  condition  of  equal  pressures  at  the  separation  points, 
only  two  parameters  of  the  physical  plane  will  be  needed  to  determine 
the  five  parameters  in  the  untransformed  plane.  These  two  parameters 
are  specified  by  the  locations  of  the  separation  points  on  the  rotating 
cylinder. 

The  first  conditions  imposed  are  designed  to  insure  that  stream¬ 
line  curvatures  and  pressures  at  the  separation  points  are  physically 


27,  L.  C.  Woods,  " Two-Dimensional  Flow  of  a  Compressible  Fluid  past 
given  Curved  Obstacles  with  Infinite  Wakes 3 "  Proc,  Ron,  So  a, 
Lond.,  Vol.  A2273  1955t  pp.  367-387. 
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meaningful  and  are  consistent  with  boundary- layer  theory.  The  curva¬ 
tures  of  the  free  streamlines  at  the  separation  points  are  required  to 
be  finite.  If  the  curvature  of  the  separation  streamline  is  negative 
infinite,  the  separation  streamline  will  be  concave  downwards  and  will 
cut  into  the  rear  of  the  cylinder.  The  negative  infinite  curvature  of 
the  separation  streamline  should,  therefore,  be  forbidden  since  this  is 
physically  impossible.  If  the  separation  streamline  has  positive 
infinite  curvature  at  .V ,  a  very  large  adverse  pressure  gradient  will 

occur  upstream  of  separation  becoming  infinite  at  Sj.  Clearly  this  is 

not  consistent  with  the  results  of  boundary  layer  theory  since  boundary- 
layer  separation  would  then  occur  before  the  assumed  separation  point. 
Thus,  only  the  case  of  finite  streamline  curvature  will  be  considered. 

The  other  condition  of  equal  pressures  at  the  separation  points 
results  from  vorticity  transport  considerations.  In  wing-airfoil 
theory,  no  net  vorticity  is  shed  from  the  wing  in  steady  flow;  that  is, 
the  average  flux  of  vorticity  out  of  any  fixed  circuit  around  the  air¬ 
foil  is  zero.  For  the  rotating  cylinder,  the  average  flux  of  vorticity 
out  of  any  fixed  circuit  around  the  cylinder  should  also  be  zero.  Now, 
if  no  vorticity  were  transported  into  the  wake  from  the  rear  part  of 
the  rotating  cylinder  immersed  in  the  wake,  then  equal  amounts  of 
vorticity  (but  with  opposite  signs)  should  be  transported  into  the  wake 
at  the  separation  points  from  the  upper  and  lower  parts  of  the  cylinder. 
The  rate  of  vorticity  transport  downstream  at  a  separation  point  is 


j  u  0 dy  •  (ues  -  uw>  ' 1  •  (ll,) 

where  the  subscript  s  denotes  conditions  at  separation.  This  is  the 
result  obtained  considering  either  the  lower  or  upper  side  of  the 
cylinder.  Since  vorticity  transport  must  then  be  the  same  at  both 
separation  points,  this  implies  that  the  magnitudes  of  ug  at  both 

separation  points  will  be  equal.  From  Bernoulli’s  equation,  the  pressure 
at  both  separation  points  will  be  the  same. 

The  original  assumption  of  equal  pressures  at  the  separation  points 
was  attributed  to  Howarth2®  by  Piercy,  Preston,  and  Whitehead.  Although 
this  assumption  is  also  used  in  the  present  work,  Piercy,  Preston,  and 
Whitehead2*  took  into  account  the  possibility  that  vorticity  is  gener- 


29.  L.  Houarth ,  "The  Theoretical  Determination  of  the  Lift  Coefficient 
for  a  Thin  Elliptic  Cylinder Proceedings  of  the  Royal  Society , 
Vol.  A-149,  1935;  pp.  558-586. 
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ated  by  the  part  of  the  cylinder  that  is  immersed  in  the  wake.  He 
obtained  good  results  for  an  ellipse  at  a  nonzero  angle  of  attack.  To 
take  this  castoff  vorticity  into  account  would  complicate  the  problam 
to  be  addressed  here  although  it  would  be  a  suitable  refinement  to  be 
considered  in  subsequent  investigations. 

C.  Numerical  Method  of  Obtaining  the  Model 

To  obtain  the  values  of  the  five  basic  parameters  appearing  in 
equation  (8)  in  the  ;  plane  -  given  the  two  separation  points  in  the 
physical  plane  -  a  system  of  five  nonlinear  equations  must  be  solved. 

One  might  attempt  to  solve  these  equations  analytically  as  was  done 
for  Bluston  and  Paulson's  relatively  simple  model.  For  their  model, 
ug  finally  reduced  to  a  trigonometric  expression  that  involved  only 

the  untransformed-plane  angles  $  and  a.  An  unsuccessful  attempt  was 
made  to  solve  the  present  more  complicated  set  of  equations  analytically; 
therefore,  it  was  decided  to  solve  for  the  five  parameters  systemat¬ 
ically  using  numerical  techniques. 

The  system  of  five  equations  is  obtained  from  requirements  im¬ 
posed  by  separation  and  the  previously  discussed  curvature  and  pres¬ 
sure  conditions  at  separation.  The  first  requirement  is  derived  from 
the  nature  of  the  critical  points,  Sj  and  S2;  these  are  the  stagnation 

points  in  the  untransformed  plane  caused  by  the  upstream  flow  from  the 
sources  placed  on  the  circle.  Thus,  the  untransformed  velocity  q  at 
the  critical  points  is  zero. 

From  the  expression  for  q,  equation  (8),  the  following  two  equations 
are  then  obtained  at  the  critical  points,  4>*±  a-e,  respectively: 

2  sin  (a  -e)  -  cot  (°  ^  '£  )  +  Q2  cot  (a  y  1  E)  J  +  uc  =  0.  (17) 

2  sin  (a  +  e)  -  £  Q1  cot  (°  )  ♦  Q2  cot  (°.  v,—  -)j-  u£  =  0.  (18) 

The  next  requirement  imposed  was  that  the  magnitude  of  the  veloci¬ 
ties  should  be  equal  at  separation  in  the  physical  plane  as  required 
by  Howarth's  assumption;  i.e: 

ue(xsl>  =  ’  (19) 

where  the  subscript  symbols  si  and  s2  indicate  the  upper  and  lower 
separation  points  respectively.  The  negative  sign  occurs  in  this 
equation  because  x  increases  on  the  upper  side  of  the  cylinder  and 
decreases  on  the  lower  side  of  the  cylinder.  Now  from  the  expression 
for  ug  in  equation  (13),  ug  is  indeterminate  at  the  separation  points 

since  q  and  N  are  both  zero  there.  This  indeterminancy  could  have 
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been  removed  if  the  problem  had  been  solved  analytically.  Since  the 
computer  was  available,  an  approximation  to  the  requirement  of 
equation  (19)  was  obtained  using  numerical  computation  techniques.  For 
angular  positions  very  near  to  the  critical  points,  the  absolute  values 
of  the  velocities  in  the  physical  plane  were  set  equal  to  each  other  at 
the  following  two  points: 

±oij  =  (a  +  Aa)  ,  (20) 

where  are  the  values  of  0  very  near  to  i  a  and  A a  is  a  small 
incremental  angle.  Now  from  equation  (15)  it  is  seen  that 

N  (otj-e)  =  N  (-a^-e) 

where  e  is  a  parameter  for  N.  Hence,  a  condition  almost  equivalent  to 
equation  (19)  can  be  imposed: 

q  (ctj  -  e)  *  -  q  (-a^  -  e)  • 

Thus,  from  equation  (8),  the  following  equation  will  result: 

r  Oj-Y-e  o^+Y+e  n 

2  [  sin  (Oj-e)  -  sin  (01^+  e)  ]  -  cov  ( — j— ■ -)  -  cot  (  j——)  J 

[a.-6-e  a.+6+e  "| 

cct  (-i^ - •)  -  cot  (—3 - •)  J  ♦  2  uc  =  0  .  (25) 

The  smoothness  or  finite  curvature  requirement  is  used  to  obtain 
two  more  equations.  A  result  of  the  smoothness  requirement  is  that  the 
value  of  du  /dx  is  finite  at  the  separation  point.  Now  from  the  chain 


rule: 

du  du  . 

ar-ar'S  •  (24) 

But  x  as  a  function  of  $  at  the  separation  points  (S^,  S^)  are  extremums 

since  the  circle  E  in  the  untransformed  plane  (Figure  2)  is  mapped  onto 
the  curved  slit  SjAS^S,  in  the  physical  plane  (Figure  3).  Therefore 


S1*S2 


and  thus 


srs2 
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is  a  necessary  condition  that  du  /dx  be  finite  at  the  critical  points. 

c 

But  N($),  which  is  zero  at  the  critical  points,  occurs  in  the  de¬ 
nominator  in  the  expression  for  du  /d$  as  a  squared  term.  So  then  in 

a  similar  manner  as  was  done  for  the  equal  pressure  condition,  use  points 
at  angle  Ao  away  from  the  points  S^,  S2: 


due(ar£)  du^-aj-e) 

d<fi  d$ 


(26) 


The  smoothness  condition  at  <p  =  -e  becomes 


m  2 

(1-2  cos  a  cos  a1  +  cos  a) 


2  (cos  a  -  cos  ap 

r  2  ai 

+  h  Qj  esc  (— 

2  - 

sin  a.  (cos  a-1)  (  f 

-  0  \  2  sin  (eye)  -  I  Qj 


2  cos  (a1  -  e) 


-e-y  ,  a  -e-6  1 

^ - ■)  +  Q2  esc  (-ij - •) 


a  -e-y 
cot  (-^ - ■) 


2  (cos  a  -  cos  a^)  ( 

+  q2  cot  (-±2 - ■)]  ♦  0  . 

The  smoothness  condition  at  <p  =  -a^  -e  is  given  as 


2  - 

(1-2  cos  a1  +  cos  a)  (  V  2  ai+E+T 

- - -  2  cos  (aye)  +  \  Qysc  ( - -■) 

2  (cos  a  -  cos  a.)  I  L  2 


(27) 


+  Q2csc  (■ 


sin  a^  (cos  a-1) 

2  (cos  a  -  cos  a^) 


2  j  -2  sin  (aye) 


aye+6  l 

+  Qr  cot  (-^ - ■)  +  Q2  cot  (~ - ■)  +  u£  j  =  0 


aye+y 


(28) 


Equations  (17),  (18),  (23),  (27),  and  (28)  constitute  a  nonlinear 
system  of  equations  to  be  solved  for  Q^,  Q2,  y,  5,  and  U£.  These  values 

are  determined  by  only  two  parameters  in  the  untransformed  plane,  a  and 
e.  These  two  parameters  are  in  turn  determined  by  x  .  and  xg2  as 
given  by  equation  (33). 
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Equations  (17),  (18),  (23),  (27),  and  (28)  are  solved  using  the  iterative 
Newton-Raphson  method  and  the  parameters  are  said  to  be  found  when  the 
biggest  change  in  any  of  the  parameters  is  less  than  2.10"^  from  one 
iteration  to  the  next. 

To  approximate  the  values  of  the  parameters  that  would  be  obtained 
if  the  pressure  and  curvature  conditions  could  be  applied  at  the  sepa¬ 
ration  points,  the  following  procedure  is  followed.  The  value  of  Ao  « 
-0.08°  is  first  used  and  the  corresponding  parameters  are  computed.  The 
parameters  for  Ao  =  0.08°  are  next  found  and  these  two  sets  of  parameters 
are  used  to  find  the  estimated  parameters  for  Ao  =  0  by  linear  inter¬ 
polation.  A  question  might  occur  as  to  how  accurate  this  approximation 
mi ght  be,  especially  if  the  parameters  changed  in  a  nonlinear  manner 
and  rapidly  with  Ao.  Figures  4  and  5  show  the  variation  of  and  y 

versus  Ao  and  both  are  approximately  linear  with  Ao  which  is  very  en¬ 
couraging  in  light  of  the  linear  interpolation  used  in  the  program. 

The  circulation  was  also  found  to  have  similar  linear  behavior. 

To  numerically  calculate  the  solution,  some  relationships  and  quanti¬ 
ties  need  to  be  found, such  as  the  radius  of  the  cylinder  in  the  z 
plane.  Using  the  mapping  function  M(0  described  by  equation  (10),  it 
is  found  that  the  position  of  in  the  z  plane  is  given  by  2i.  Drawing 

a  line  from  the  center  of  the  circular  arc  C  in  the  z  plane  to  Sj,  the 

radius  is  seen  to  be 


a  =  2  esc  m  ,  (29) 

si 

where  --as  mentioned  before--the  subscript  si  denotes  evaluation  at  the 
upper  separation  point  S^.  The  length  from  the  origin  in  the  z  plane 

along  the  real  line  to  the  arc  C  is  given  by 

M  (-esc  a)  =  -  2  cot  a  .  (30) 

Thus,  drawing  a  straight  line  from  point  A  to  Sj ,  it  is  seen  from 
equation  (30)  that 

a  =  /2  .  (31) 

Using  equation  (31),  a  is  then  found  to  be,  from  equation  (29): 

a  =  2  esc  2  a  .  (32) 

Knowing  and  x^,  the  values  a  and  e,  which  determine  the  five 

parameters  in  the  plane  r,,  can  be  found  using  Figure  3  and  equation 
(31): 

E  =  (xsl  +  Xs2)/2  ’  (33) 
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Act  (degrees) 

Figure  5.  y  vs.  a a 
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(33) 


«  B  (»  -  Xsl+E)  /2 

Since  the  quantities  in  the  plane  z  will  be  calculated  in  terms 
of  the  angular  position  $  in  the  ?  plane,  it  is  necessary  to  determine 
a  relationship  between  tf>  and  x.  To  do  this,  consider  the  directed 
vector  t  from  the  center  of  the  circular  arc  C  in  the  physical  planes 
to  a  point  on  the  arc  C  in  the  physical  planes.  Using  the  mapping 
function  given  by  equation  (10) ,  the  expression  for  ?  is  given  by  the 
following: 

a  =  M  I  esc  a  exp  (io)]  -M(0)  .  (34) 

This  expression  can  be  manipulated  to  a  form  such  that  the  imaginary 
part  of  t  can  be  taken.  The  imaginary  part  of  a  is  just  a  sin  m.  Using 
the  imaginary  part  of  the  expression  for  a  from  equation  (34)  and 
applying  trigomonetric  identities,  the  following  relationship  is  obtained: 


sin  a  (1-cos  a  cos  o) 
sin  m  =  . -  . f 

(sec  a  +  1/sec  a)  -  cos  o  .  (35) 


Expressing  equation  (35)  in  terms  of  $  and  x,  the  following  equation  is 
obtained: 


sin(x-e) 


lt£J . Ik— a  .c.°4i»l0l.  (361 

[(sec  a  +  1/sec  a)  /2J-  cos  (*+e)  1  ' 


The  velocity  and  derivative  of  the  velocity  with  respect  to  angle 
in  the  physical  plane  are  calculated  as  these  quantities  appear  in  the 
integral  boundary- layer  equations.  A  subprogram  was  constructed  to 
obtain  these  quantities  found  in  terms  of  the  untransformed  variable 
(if)  and  parameters  in  the  plane  ?.  This  subprogram  is  called  at  each 
step  of  the  integration  of  the  boundary- layer  equations.  First  the 
relationship  between  x  and  $  is  found  using  equation  (36) .  The  value 
of  $  is  found  using  the  iterative  method  of  Newton.  Successive  values 
of  <|>  found  during  iteration  must  be  within  10"®  before  the  method  is 
declared  to  have  converged.  The  values  of  u  and  du  /dx  can  then  be 

found  from  equation  (8),  its  derivative  with  respect  to  $,  equation 
(15),  its  derivative  with  respect  to  4> >  and  substituting  the  results 

into  equations  (13)  and  its  derivative  with  respect  to  For  more 
details,  the  listing  of  the  computer  program  is  given  in  Appendix  A. 


D.  Method  of  Obtaining  Agreement  Between  the  Flow  Model  and  the 
Boundary-Layer  Calculation  Results 

For  bodies  with  fixed  surfaces,  boundary- layer  calculations  can 
approximately  predict  the  point  of  separation  given  the  external 
velocity  distribution.  In  this  investigation  of  the  rotating  cylinder, 


31 


the  final  inviscid-flow  model's  separation  points  are  made  to  agree 
with  the  location  of  the  separation  points  found  by  numerical  calcu¬ 
lation  of  the  boundary  layer.  This  flow-model  can  thus  be  made  con¬ 
sistent  with  boundary- layer  theory  using  the  integral  boundary- layer 
technique  developed  here. 

The  procedure  for  obtaining  consistency  between  the  flow-model 
and  boundary  layer  theory  is  quite  simple  in  principle.  The  boundary 
layer  equations  for  flow  over  the  rotating  circular  cylinder  are 
solved  with  assumed  separation  points  for  the  inviscid-flow  model 
further  back  on  the  cylinder  than  they  would  be  anticipated  to  actually 
occur.  Separation  points  are  then  obtained  from  the  results  of  the 
boundary- layer  calculation  and  a  revised  inviscid  flow  model  is  then 
constructed  using  these  new  separation  points.  Again  new  boundary- 
layer  separation  points  are  found  and  the  flow  model  is  again  modi¬ 
fied  using  these  new  separation  points.  This  process  is  repeated 
until  the  largest  absolute  value  of  the  difference  in  separation 
location  for  the  last  two  particular  flow  models  is  less  than  a  value 

called  E  .  The  value  of  E  used  in  this  work  was  E  =0.46  degrees, 
c  c  c 


III.  FLOW-MODEL  RESULTS  AND  DISCUSSION 

The  external  velocity  distributions  consistent  with  the  boundary- 
layer  calculations  and  the  corresponding  lift  and  drag  coefficients  are 
of  chief  interest  for  comparison  with  experiment.  The  range  of  uw 

considered  here  (0  <  u  <  0.32)  is  chiefly  limited  by  stability  of  the 

w  1 

boundary- layer  calculation  method  used  although  Swanson  indicates  that 
turbulence  always  appears  before  separation  on  the  lower  side  of  the 
cylinder  when  uw  >  0.5  for  the  range  of  Reynolds  numbers  investigated. 

He  investigated  the  Magnus  effect  for  3.58-10^  <  Re^  <  5.01°10 


A.  External-Velocity  Distributions 

Figure  6  illustrates  the  first  trial  external -velocity  distribution 
up  to  the  point  of  separation  and  compares  it  with  the  converged  or 
final  velocity  distribution  when  uy  =  0.2.  Figures  7  and  8  present,  in 

the  adverse  pressure  gradient  region,  details  of  these  distributions  with 
additional  distributions  found  in  the  iteration  process  to  obtain  con¬ 
vergence.  The  graphs  are  representative  of  the  convergence  behavior 
using  this  iterative  method  for  all  corresponding  values  of  uw  considered 

with  one  exception.  For  uw  =  0.3,  various  initial  values  of  separation 

points  were  tried,  but  after  a  few  iterations  the  boundary  layer  calcu¬ 
lations  would  not  separate  up  to  and  including  the  breakaway  point  of 
the  inviscid-flow  model.  It  is  felt  that  the  boundary  layer  needs  to 
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X  (DEGREES) 

Figure  7.  Velocity  Distribution  Details  -  Upper  Side 
Adverse  Pressure  Gradient  Region 
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be  calculated  more  accurately  near  separation  for  iarger  negative  ^ 
values  of  uw/u@.  Granted  adequate  time,  it  is  thought  that  consisten 

and  more  meaningful  results  could  be  achieved  by  improving  the  integral 
boundary- layer  technique. 


For  a  cylinder  rotation  rate  of  0.15,  the  locations  of  separation 
versus  the  trial  number  in  the  iteration  process  to  find  a  converged 
solution  is  shown  in  Figure  9.  As  mentioned  before,  when  the  absolute 

value  E  of  the  differences  between  the  positions  of  succeeding  values 
c 

of  separation  points  on  both  sides  of  the  cylinders  become  less  than 

0.46  degrees,  convergence  was  considered  to  be  achieved.  Bluston  and 

Paulson's  value  of  E  was  0.5°.  For  u  =  0.15,  the  value  of  E  was 

c  w  C 

decreased  to  0.16  degrees  to  observe  the  convergence  behavior  more 
closely.  In  Figure  9,  it  is  seen  that  three  more  iterations  were  re¬ 
quired  to  satisfy  the  tighter  convergence  criterion. 

The  convergence  of  the  lift  coefficient  can  also  be  studied  in 
Figure  10.  Here  it  is  seen  that  the  value  of  CL  from  circulation 

considerations  is  quite  sensitive  to  separation  location.  The  value 

of  C.  for  E  =0.16  was  about  80%  of  C.  for  E  =  0.46  degrees.  However, 
L  C  L  c 

by  integrating  the  vertical  component  of  the  pressure  coefficient  over 

the  cylinder  up  to  the  point  of  separation,  the  values  of  CL  were  found 

to  be  within  10%  of  each  other  for  the  two  values  of  Ec> 


In  choosing  a  precision  criterion,  one  must  accept  a  compromise 
between  computing  time  and  accuracy  of  the  boundary  layer  calculation 
method.  In  retrospect,  it  appears  that  by  using  an  accelerated  con¬ 
vergence  technique,  a  more  rigorous  standard  for  the  convergence  cri¬ 
terion  could  have  been  adopted.  However,  it  might  be  better  to  im¬ 
prove  the  inviscid  flow  model  in  some  other  wayj  for  instance,  one 
could  change  the  flow  model  to  take  into  account  the  cast-off  vorticity 
from  that  part  of  the  cylinder  immersed  in  the  wake  region. 

Figure  11  compares  some  experimental  velocity  distributions  with 
the  results  of  the  current  analysis  for  the  nonrotating  cylinder.  The 
curve  of  Fage  and  Falknor29  was  obtained  from  Bluston  and  Paulson's 
paper  while  the  curve  of  Petrie  and  Simpson30  was  taken  directly  from 


29.  A.  Fage  and  V.  M.  Falkner ,  "The  Flow  Around  a  Circular  Cylinder, " 
ARC  E&M  No.  1369,  1931. 

30.  A.  M.  Petrie  and  H.  C.  Simpson,  "An  Experimental  Study  of  the 
Sensitivity  to  Freestream  Turbulence  of  Heat  Transfer  in  Wakes 
of  Cylinders  in  Crossflow,"  Int.  J.  of  Heat  and  Mass  Transfer, 
Vol.  15,  1972,  pp.  1497-1513. 
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Figure  9.  Location  of  Boundary-Layer  Separation  vs.  Trial  Number 
(uw  -  0.15) 
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THEORETICAL  RESULTS 
FAGE  »  FALKNER  (R«d  - 1 05 ) 

PETRIE  8  SIMPSON  (  Red  s  I  -  4x10*  ) 


their  paper.  Fage  and  Falkner's  curve  gives  the  better  agreement  with  the 
theoretical  curve,  but  is  for  a  higher  Reynolds  number  than  for  the  other 
curves  with  the  amount  of  turbulence  in  the  fveestream  unknown.  The  curve 
of  Petrie  and  Simpson  was  obtained  with  very  little  mainstream  turbulence; 
by  introducing  a  large  amount  of  turbulence  into  the  mainstream,  they 
obtained  data  similar  to  the  Fage  and  Faulkner's  results.  The  theo¬ 
retical  curve  approximates  the  experimental  curves  better  than  the 
simple  sine  curves  or  series  curves.  Nevertheless,  some  sort  of  re¬ 
finement  of  the  inviscid-flow  model  might  better  represent  the  ob¬ 
served  pressure  distributions  in  the  adverse  pressure  gradient  region. 

A  check  of  Bluston  and  Paulson's  velocity  distribution  reveals  it 
to  be  nearly  identical  to  the  theoretical  distribution  found  in  the 
present  investigation.  This  agreement  with  Bluston  and  Paulson  further 
increases  confidence  in  the  procedure  of  finding  the  five  parameters 
(Qj.  Q,»  Y >  <5,  u£)  in  the  untransformed  plane  by  computer. 

It  has  been  seen  that  converged  separation  points  can  be  found  by 
this  iterative  process  given  certain  initial  trial  separation  points. 

The  questioii  naturally  arises  as  to  how  the  final  separation  points 
are  affected  by  the  particular  choice  of  initial  separation  points. 

To  answer  this  question,  different  sets  of  initial  separation  point 
values  were  used  in  the  computer  program.  The  results  are  tabulated 
as  follow: 


TABLE  I. 

INITIAL  AND 

FINAL  VALUES  OF 

SEPARATION 

Initial 

Values 

Final 

Values 

Upper 

Side 

Lower 

Side 

Upper 

Side 

Lower 

Side 

u 

w 

Xsl 

Xs2 

x  i 
si 

Xs2 

0.2 

95.1° 

-81.9° 

87.50° 

-70.61° 

0.2 

92  8° 

-83.1° 

87.57° 

-70.82° 

0.1 

94.5° 

-87.6° 

85.77° 

-77.30° 

0.1 

95.0° 

-81 .9° 

85.86° 

-77.28° 

The  final  values  for  the  separation  points  vary  at  most  0.2°  with 
these  initial  choices.  As  will  be  seen  later,  the  variation  in  the 
lift  and  drag  coefficients  will  also  vary  a  small  amount  with  initial 
choice  of  the  breakaway  points.  This  dependence  of  the  final  values 
on  the  initial  values  could  probably  be  decreased  further  by  using  a 
more  rigorous  standard  for  convergence. 
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Figure  12  shows  absolute  values  of  u  plotted  against  the  absolute 

value  of  x  for  some  different  values  of  uw.  The  separation  points  can 

be  seen  to  be  delayed  on  the  upper  side  and  advanced  on  the  lower  side 
in  agreement  with  observation. 

The  inviscid  flow  fields  external  to  the  cylinder  were  also  in¬ 
vestigated.  Figure  13  shows  the  streamlines  in  the  untransformed  plane 
for  the  nonrotating  cylinder  and  Figure  14  exhibits  the  streamlines  in 
the  corresponding  physical  plane.  Figure  15  and  16  gives  the  pattern 
of  streamlines  in  the  untransformed  and  physical  plane  respectively 

for  u  =  0.2.  The  zero  stream  function  lines  in  the  untransformed 
w 

plane  are  seen  to  leave  perpendicular  to  the  cylinder  whereas  in 
the  transformed  plane  these  lines  leave  the  circle  tangentially  as 
required  by  the  transformation.  It  may  be  observed  that  the  stream¬ 
lines  that  form  the  boundary  of  the  wake  in  the  physical  plane  are 
very  nearly  parallel;  the  main  effect  of  the  low  rotation  rate 
is  to  produce  a  slight  asymmetry  in  the  streamlines  near  the  separa¬ 
tion  points.  The  streamline  pattern  in  the  untransformed  plane  show 
greater  departures  from  a  symmetrical  pattern  than  for  the  physical 
plane  case. 

The  stagnation  points  shift  to  negative  x  values  as  uw  is  in¬ 
creased.  Below  is  a  table  giving  these  stagnation  points  for  dif¬ 
ferent  values  of  u  . 

w 


*The  curves  do  not  extend  to  x  =  0  since  values  of  u  and  du  /dx  are 

e  e 

printed  out  only  for  that  region  being  numerically  calculated  with 
the  integral  boundary -layer  equations.  The  limitations  of  the 
computer  code  restrict  integration  of  the  boundary -layer  equations 
on  the  lower  side  to  the  region  where  u yu&  >_  -0. 5.  Thus ,  the  value 

of  |x|  at  which  integration  can  start  becomes  larger  with  increasing 
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Figure  13.  Streamline  Pattern  Around  Cylinder  in  Untransformed 
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Fitiur*  14,  streamline  PaLtum  Around  Cyl Under  In  Physical  Plane 


Hgui-t  16.  Streamline  Pattern  ArouruJ  Cylinder  iri  Pny&lcit  Plane 


TABLE  II.  STAGNATION  POINT  ON  CYLINDER  FOR  DIFFERENT 
RATES  OF  ROTATION 

uw  Stagnation  Point  Location  (Degrees) 


0.0 

0.0 

0.1 

-0.53 

0.2 

t 

o 

to 

t— * 

0.25 

-1.19 

0.32 

-1.76 

This  displacement  direction  of  the  stagnation  point  also  occurs  in 

experiments.  However,  no  detailed  comparison  with  experiment  can  be 

made  for  these  small  rotation  rates  partly  because  of  the  difficulty 

in  measuring  small  changes  in  the  stagnation  point  location.  Below 

u  =0.25,  the  stagnation  point  is  displaced  linearly  with  u  but  the 
w  w 

displacement  for  uw  =  0.32  is  larger  than  the  linear  relationship  would 

predict.  The  velocity  distribution  curves  uw  =  0.2  and  uw  =  0.25  do 

not  follow  the  trends  expected  and  noted  in  Figure  12.  Comparing  the 
upper-side  velocity  distributions  in  Figure  17  reveals  that  the  maximum 
value  of  velocity  for  uw  =  0.15  is  actually  higher  than  it  is  for 

uw  =  0.25.  The  upper  side  displacement  of  separation  increases  only  a 

minimal  amount  with  the  increase  in  uw<  However,  the  lower  side  curve 

follows  the  trend  noted  for  the  lower  values  of  wall  velocity.  The 
reasons  for  this  behavior  are  not  known. 

For  laminar  boundary  layers,  the  positions  of  separation  should  be 
a  monotonic  function  of  uw  according  to  experimental  indications  and 

also  as  expected  from  boundary  layer  theory.  Figure  18  shows  the  com¬ 
puter  results  for  the  absolute  values  of  the  separation  points  versus 
uw»  It  appears  that  for  values  of  uw  <  0.15,  the  displacements  of  the 

separation  points  from  the  nonrotating  case  are  almost  linear.  However, 

the  absolute  value  of  the  slope  for  the  line  approximating  the  computer 

results  is  greater  for  the  lower  side  than  it  is  for  the  upper  side  of 

the  cylinder.  In  prior  work,^  where  sine- function  external-velocity 

distributions  were  used  with  this  integral  method,  a  linear  variation 

with  u  was  also  observed  but  extending  up  to  values  of  u  =  0.4.  The 
w  w 

slope  of  this  line  is  somewhat  less  than  for  the  present  bound-vortex 

source-wake  model. 
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|X|  (DEGREES) 

Figure  17.  |ue|vs.|x|—  uw  =  0.25  and  uw  =  0.15 
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B.  Lift  and  Drag  Coefficients 


Approaches  representing  two  different  views  can  be  utilized  to  ob¬ 
tain  the  lift  coefficient.  The  first  approach  stems  from  the  consider¬ 
ation  that  the  rotating  cylinder,  through  the  action  of  the  boundary 
layer,  has  induced  in  the  inviscid  flow  a  certain  circulation  strength 
around  the  cylinder.  This  inviscid  flow  model  takes  this  induced  cir¬ 
culation  into  account  by  introducing  a  bound  vortex  in  the  physical 
plane  or  more  precisely,  a  point  vortex  of  strength  r  at  the  center  of 
the  circle  in  the  untransformed  plane.  This  bound  vortex  gives  a  lift 
force  according  to  the  Zhukovskii  theorem.  In  the  other  approach,  one 
is  concerned  with  providing  a  flow  model  with  velocity  distributions 
approximating  the  experimental  distributions.  The  assumption  of 
Howarth26  is  then  made  in  which  constant  pressure  is  assumed  over  that 
part  of  the  cylinder  immersed  in  the  wake.  The  lift  coefficient  can 
thus  be  calculated  knowing  the  velocity  distribution.  No  attempt  will 
be  made  here  to  reconcile  these  two  approaches  or  choose  between  them. 
Nevertheless,  it  might  be  pointed  out  that  the  latter  approach  is  the 
same  approach  used  in  obtaining  the  drag  coefficients.  Calculations 
were  made  for  both  approaches:  the  bound- vortex  approach  and  the  in¬ 
tegrated-pressure  approach. 

The  lift  coefficient  for  the  bound- vortex  approach  can  be  found 
by  first  considering  that  the  circulation  strength  and  hence  the 
lifting  force  of  the  point  vortex  can  be  shown  to  be  invariant  under 
the  transformation  between  the  two  complex  planes.  The  lift  co¬ 
efficient  can  then  be  found  to  be 

C,  =  2ir  u  cos  a 
L  c 

In  the  second  approach,  the  lift  coefficient  is  obtained  by  first 
integrating  the  vertical  component  of  the  local  force  coefficient  on 
the  cylinder  as  calculation  of  the  boundary- layer  on  the  cylinder  pro¬ 
ceeds.  The  vertical  component  of  the  force  coefficient  on  the  cylin¬ 
der  surface  in  the  wake  is  then  calculated  assuming  the  pressure  con¬ 
stant  in  the  region.  Finally,  these  contributions  are  used  to  ob¬ 
tain  the  total  lift  force  coefficient  on  the  cylinder. 

Using  these  two  approaches,  the  lift  coefficient  as  a  function 
of  rotation  rate  (Ow/uq)  is  shown  in  Figure  19.  These  lift  coef¬ 
ficient  values  are  compared  with  the  experimental  results  obtained  by 
Swanson  for  Reynolds  numbers  from  3.58x10^  to  9.9x10  .  The  integrated 
pressure  approach  gives  the  better  agreement  for  most  of  the  range  of 

u  and  the  lift  coefficient  is  linear  with  u  up  to  and  including 
w  ™ 

u  =  0.15.  For  u  =  0.2,  two  computer  runs  were  made  with  different 
w  w 

initial  values  of  separation  points;  the  variation  in  the  final  results 
are  also  shown  in  Figure  19.  The  variation  of  the  lift  coefficient 
with  variation  of  initial  v  iues  could  probably  be  lessened  by  making 
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Figure  19.  C,  vs.  uw  -  Comparison  with  Experiment 


the  convergence  criterion  for  this  iteration  procedure  more  severe. 


To  obtain  the  drag  coefficient,  the  horizontal  component  of  the 
local  force  coefficient  on  the  cylinder  is  integrated  as  the  boundary 
layer  quantities  are  calculated;  the  contribution  to  the  drag  co¬ 
efficient  in  the  wake  region  is  found  assuming  that  it  is  a  region  of 
constant  pressure.  The  drag  coefficient  obtained  by  computer  is  shown 
together  with  experimental  results  in  Figure  20.  The  values  obtained 
from  the  model  are  comparable  to  values  found  experimentally,  although 
the  general  trends  do  not  agree  very  well.  Similarly,  as  for  the  lift 
coefficient  comparison,  the  value  at  uw  =  0.25  is  in  the  worst  agree¬ 
ment  with  the  experimental  value. 

The  values  of  the  obtained  untransformed-plane  parameters 
corresponding  to  the  lift  and  dug  force  results  are  summarized  in 
the  following  table. 


TABLE  HI.  PARAMETERS  OBTAINED  FOR  THE  UNTRANSFORMED  PLANE 
(Angles  are  in  radians  except  that  values  in 
parentheses  are  expressed  in  degrees) 


u 

w 

<*1 

q2 

Y 

8 

u 

c 

a 

f 

0.0 

.28918 

.28918 

.3322 

-.3322 

0.0 

.8551 

0.0 

(48.994) 

(0.0) 

0.1 

. 18564 

.40316 

.32267 

-.34984 

.01125 

.85898 

.07492 

(49.216) 

(4.2926) 

0.15 

.13037 

.48559 

.3324 

-.35289 

.0132 

.87119 

.11381 

(49.916) 

(6.5208) 

0.2 

.082203 

.55964 

.3552 

-.3545 

.016 

.88089 

.14748 

(50.471) 

(8.45) 

0.25 

.055223 

.  64684 

.41081 

-.34946 

.02198 

.89503 

.18146 

(51.281) 

(10.397) 

0.52 

.08099 

.49354 

.3164 

-.4011 

.0484 

.84293 

.1744 

(48.296) 

(9.9924) 

Here  it  is  seen  that  the  angles  y  and  6  are  relatively  insensitive  to 
separation  point  locations.  The  reasons  for  the  insensitivity  arc  not 
known. 

This  flow  model  could  probably  be  improved  by  taking  into  account 
the  vorticity  shed  from  that  part  of  the  cylinder  immersed  in  the  wake 
This  is  deduced  from  Swanson's  observation  that  for  uw  =  0.2  and 
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Figure  20.  Cn  vs.  u,.,  -  Comparison  with  Experiment 


4 

Re^  =  4x10  ,  the  pressure  at  the  upper  separation  point  is  lower  than 

it  is  for  the  separation  point  where  the  wall  is  moving  upstream. 
According  to  equation  (16)  and  Bernoulli's  equation,  this  implies  that 
the  transport  of  vorticity  into  the  wake  at  the  upper  separation  point 
is  greater  than  for  the  vorticity  of  opposite  sign  transported  into 
the  wake  at  the  lower  separation  point.  Thus,  vorticity  of  the  same 
sign  as  for  the  lower  cylinder  surface  must  be  cast  off  from  the 
cylinder  surface  immersed  in  the  wake.  This  feiture  might  be  taken  in¬ 
to  account  by  setting  the  velocity  magnitudes  of  the  two  free  stream¬ 
lines  equal  to  each  other  where  they  intersect  a  vertical  line  perhaps 
three  or  four  diameters  downstream  of  the  cylinder.  For  an  ellipse 
at  angle  of  attack,  Piercy,  Preston  and  Whitehead  took  the  cast-off 
vorticity  in  the  wake  region  into  account  by  imposing  the  equal  veloc¬ 
ity  condition  along  two  vertical  lines  cutting  the  wake  a  distance  c 
and  2c  from  the  ellipse.  Here  c  is  the  length  of  the  long  axis  of  the 
cylinder.  Differences  in  lift  between  calculation  and  observation 
were  at  most  15%,  with  angle  of  "stall"  for  the  ellipse  being  predicted 
with  good  accuracy. 


IV.  SUMMARY  AND  CONCLUSIONS 

An  inviscid  flow-with-wake  model  for  a  spinning  cylinder  in  cross- 
flow  was  developed  that  is  consistent  with  boundary- layer  theory.  This 
was  accomplished  by  first  considering  flow  around  a  cylinder  in  an  un¬ 
transformed  plane  with  two  sources  on  the  cylinder,  their  image  sinks 
at  the  center  and  a  point  vortex  in  the  center  of  the  circle.  This 
resultant  flow  was  transformed  to  the  physical  plane  by  a  Zhukovskii 
transformation;  this  results  in  the  separation  streamlines  at  the 
separation  points  being  tangent  to  the  cylinder  at  the  critical  points 
of  the  transformation. 

The  condition  that  the  pressure  gradients  are  finite  at  the  sepa¬ 
ration  points  was  imposed  on  this  elementary  inviscid  flow-model.  The 
imposed  pressure-gradient  condition  produces  a  more  realistic  inviscid 
flow  near  separation.  This  condition  also  makes  possible  the  develop¬ 
ment  of  a  self-consistent  model;  that  is,  the  inviscid  flow  model  takes 
the  boundary- layer  into  account  and  is  consistent  with  boundary- layer 
theory.  With  finite  pressure  gradients  at  separation,  the  number  of 
parameters  needed  to  describe  the  flow  was  reduced  from  five  to  three. 

Another  condition  imposed  upon  the  inviscid  flow  model,  equal 
pressures  at  the  separation  points,  results  from  vorticity  transport 
considerations.  The  vorticity  transport  into  the  wake  at  the  separa¬ 
tion  points  was  assumed  to  be  equal  with  no  vorticity  cast  off  from 
the  rear  of  the  cylinder.  After  applying  this  condition,  only  two 
parameters  of  the  flow  in  the  physical  plane  were  needed  for  the 
description  of  this  flow.  The  two  parameters  used  were  the  location 
of  the  separation  points  on  the  cylinder.  Using  these  parameters,  a 
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system  of  five  equations,  which  describes  these  imposed  conditions, 
was  solved  for  the  five  basic  parameters  in  the  untrans formed  plane. 
These  equations  were  solved  numerically  by  computer. 

An  iteration  procedure  was  used  to  obtain  separation  points 
located  in  agreement  with  boundary- layer  theory.  The  boundary- layer 
was  calculated  using  the  integral  technique  and  new  separation  points 
were  obtained  using  a  particular  inviscic-flow  model  with  assumed 
separation  points.  These  new  separation  points  found  by  the  integral 
method  then  become  the  separation  points  for  a  new  assumed  inviscid- 
f low  field  and  again  the  boundary  layer  is  calculated  to  find  a  new 
separation  point.  The  iteration  procedure  was  stopped  when  the 
difference  between  two  successive  separation  points  was  less  than  0.46 
degrees  for  both  sides. 

Comparison  of  the  converged  velocity  distribution  for  the  non¬ 
rotating  case  with  that  of  the  observed  velocity  distributions  shows 
fair  agreement.  Comparison  between  the  calculated  velocity  dis¬ 
tributions  for  different  values  of  u^  showed  that  separation  is  de¬ 
layed  and  the  maximum  value  of  ug  increases  with  uw  for  the  upper  or 

downstream-moving  wall  side  of  the  cylinder  with  the  exception  of  the 
velocity  distributions  for  uw  =  0.2  and  u^  =  0.25.  For  the  lower  side 

of  the  cylinder,  or  in  other  words,  the  side  with  the  upstream  moving 
wall,  separation  advances  toward  the  front  of  the  cylinder  with  in¬ 
creasing  uw  without  exception.  The  stagnation  point  moves  downward 

on  the  cylinder  with  increasing  uw  as  is  observed  in  experiment. 

The  calculated  lift  coefficient  was  compared  with  the  observed 
lift  coefficient.  When  the  lift  coefficient  was  calculated  using  the 
integrated  value  of  the  pressure  over  the  cylinder,  a  linear  (also, 
obtained  by  experiment)  relationship  was  obtained  up  to  and  including 

u  =0.15,  and  the  lift  coefficient  increased  almost  linearly  after 
w 

uw  =  0.15.  Using  the  circulation  or  bound-vortex  approach,  the 

calculated  lift  coefficients  were  somewhat  lower  than  the  preceding 
approach  but  still  in  good  agreement  with  experiment  up  to  and 
including  =  0.15.  The  drag  coefficients  for  different  uw  were 

also  calculated  and  are  at  most  20%  lower  than  the  observed  value. 

Although  the  approach  does  have  some  success  in  predicting  the 
Magnus  force  on  a  rotating  cylinder,  refinements  and  modifications  of 
the  present  method  could  possibly  improve  the  computed  results 
significantly.  Some  refinements  and  modifications  that  could  be 
applied  to  the  present  model  are  the  following: 

(1)  tighter  convergence  criterion  and  use  of  an  accelerated 
convergence  technique, 
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(2)  modification  of  model  to  take  into  account  the  vorticity 
shed  into  the  wake  from  that  part  of  the  cylinder  which  is  immersed 

in  the  wake,  and 

(3)  use  of  a  more  complicated  distribution  of  sources. 
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APPENDIX  A 


LISTING  OF  COMPUTER  PROGRAM  TO  CALCULATE  FORCE  COEFFICIENTS 

The  listing  of  the  entire  program  is  given  here  except  for  the 
data  used  to  form  the  integral  quantity  arrays.  The  computer  program 
is  written  in  Fortran  IV.  The  available  shape  factors  are  given  in 
tabulated  form  in  Reference  24.  The  list  of  principal  variables  for 
the  program  is  given  below  for  the  main  program. 


Fortran  Implementation 
List  of  Principal  Variables 


PROGRAM  SYMBOL 


DEFINITION 


A 

AB,  ABP,  ABM 

ADELT 

AGAM 


AH,  AHP,  AHM 


similarity  solution  type,  A 

two-dimensional  arrays  of  the  pressure- 
gradient  parameter  3, 

angular  location  of  lower  source  in  the 
£  plane 

angular  location  of  upper  source  in  x, 
plane 

two-dimensional  arrays  of  the  shape 
factors  H 


AK,  AKP,  AKM 

AL,  ALP,  ALM 


ordered  one-dimensional  arrays  of  the 
shape  factors  K 

two-dimensional  arrays  of  the  shape 
factors  L 


AKL 


AKS 


AT,  ATP,  ATM,  AUS 


one-dimensional  array  of  highest  values 

of  K  corresponding  to  values  u  /u 

w  e 

represented  by  the  one-dimensional 
array  AUS 

one-dimensional  array  of  separation 
values  of  K  corresponding  to  the  values 
uw/ug  represented  by  the  one  dimensional 

array  AUS. 

two-dimensional  array  of  u  /u  values 

w  e 

used  for  interpolation  purposes. 


PROGRAM  SYMBOL 


DEFINITION 


AU 


AUM 


AUPP 


B 


one-dimensional  array  of  values  of  u  /u 

w  e 

corresponding  to  the  columns  of  the  two- 
dimensional  arrays  AH,  AT,  AL. 

one-dimensional  array  of  values  of  u  /u 

w  e 

corresponding  to  the  columns  of  the  two- 
dimensional  arrays  AHM,  ATM,  ALM 

one-dimensional  array  of  values  of  u  /u 

w  e 

corresponding  to  the  columns  of  the  two- 
dimensional  arrays  ATP,  AHP,  and  ALP. 

velocity  parameter  defined  by  B=ue/(uw-ue) 


BB 

BETA 

Bll 

BK 

BT 

BET1 ,  BET 2 


BETN1 ,  BETN2 

CIRC 

CKF 


COEFI) 


two-dimensional  array  of  values  of  B 

pressure  gradient  parameter^ 

two-dimensional  array  of  values  of  the 
shape  factors  11. 

two-dimensional  arrays  of  values  of  the 
shape  factors  K. 

two-dimensior.al  array  of  values  of  the 
shape  factors  T. 

initial  guessed  value  of  the  separation 
point  on  upper  side  and  lower  side 
respectively.  Also  used  to  store  the 
last  found  separation  values. 

most  recent  values  of  the  separation  points 
on  the  upper  and  lower  sides  respectively. 

nondimensionalized  vortex  strength. 

initial  value  of  K  to  start  the  integra¬ 
tion  of  the  boundary- layer  equations. 

drag  coefficient. 


COEEL 


lift  coefficient  calculated  from  circu¬ 
lation  considerations. 
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PROGRAM  SYMBOL 


DEFINITION 


COEFLF 


lift  coefficient  calculated  from  pressure 
considerations. 


CPC0EF  local  pressure  coefficient. 

CUI  starting  value  of  velocity  ratio,  this 

determines  the  starting  value  of  x  along 
the  cylinder. 

DELTA  the  non-dimens ionali zed  displacement 

thickness. 


DERIVX  subroutine  which  computes  the  derivaties 

and  stores  them  in  DX. 

DMAX  absolute  value  of  the  maximum  step-size 

to  be  used. 


DNXT 

DPST 

DUX 

DX 


DERIVX 

EPCON 

ER 


Fll 


step  size  to  be  used  for  next  step. 

actual  step  size  used  with  the  step 
just  completed 

derivative  of  boundary- layer  edge  velocity 
with  respect  to  position. 

one  dimensional  array  for  storing  the 
derivative  obtained  by  the  subroutine  for 
evaluating  the  derivative  . 

a  subroutine  which  computes  the  derivatives 
and  stores  them  in  DX. 

value  used  as  a  test  to  see  if  convergence 
has  been  achieved. 

amount  of  maximum  error  for  each  step. 

If  the  error  is  larger  than  ER,  the 
step  size  is  reduced. 

variable  characterizing  the  skin  friction. 


KUTMER 


variable  step  Runge-Kutta  integration 
subroutine. 


N 


number  of  equations  to  be  evaluated  in 
DERIV  subroutine. 
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PROGRAM  SYMBOL 
NBP 

NEQ 

NFIRST 

NIT 

NPAR 
PR  l  NX 
PS 

Q1 

Q2 

THATA 

THETA 1 

TERMX 

UI 

uw 

VELOC 


DEFINITION 

number  of  elements  in  the  arrays  AUS, 

AKS,  and  AKL. 

number  of  equations  to  be  solved  simultan¬ 
eously  for  the  basic  parameters  in  the 
plane.  An  input  to  the  subroutine  VELPAR. 

values  which  denotes  whether  the  upper  or 
lower  half  of  the  cylinder  is  being 
considered. 

number  of  iterations  required  to  obtain 
UI  and  DUX. 

a  number  value  equal  to  NEQ  plus  one. 
print  subroutine. 

print  step  value.  PS  determines  the 
frequency  of  reference  to  subroutine 
PRINX. 

value  that  is  one-half  of  the  upper 
dimensionless  source  strength. 

value  that  is  one-half  of  the  lower 
dimensionless  source  strength. 

the  non-dimensionalized  momentum  loss 
thickness . 

initial  value  of  the  non-dimensional 
momentum  thickness. 

termination  subroutine. 

local  external  velocity  on  cylinder. 

tangential  velocity  of  rotating  cylinder 
non-dimensionalized  by  the  free-stream 
velocity. 

subroutine  to  calculate  the  edge  velocity 
and  its  derivative  given  a  position  on 
the  cylinder. 
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PROGRAM  SYMBOL 


DEFINITION 


VELPAR 


W 

X 

XINTC 

XINTSQ 

XNEW 

XOLD 

XPV 

XTC 

xz 

XI 
X1D 
XI  INC 


subroutine  used  to  find  the  basic  param¬ 
eters  of  the  ;  plane.  It  calls  other 
subroutines  in  order  to  effect  this 
mission. 

one-dimensional  array  of  length  2N.  Not 
used  here  but  is  dimensioned  in  main 
program. 

one-dimensional  array  of  length  N,  which 
contains  the  variables  found  by  numerical 
integration. 

X 

integrated  value  of  G*  .  Needed  to  compute 
the  energy  loss  thickness  and  used  to 
obtain  the  shape  factor  values. 

2 

integrated  value  of  G"  .  Used  in  ob¬ 
taining  some  shape  factor  values. 

an  array  containing  the  most  recently 
found  values  of  the  basic  parameters 
of  the  ;  plane. 

initial  guessed  values  of  the  basic 
parameters  of  the  c  plane. 

print  variable 

one-dimensional  array  of  length  NTS. 

KUTMER  will  return  to  main  program  if 
any  element  of  XTC  changes  sign. 

position  on  circle  in  c  plane. 

position  on  cylinder  in  z  plane. 

the  value  of  XI  in  degrees. 

integration  step  size  used  to  find  force 
coefficients  before  KUTMER  subprogram  is 
used. 
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r  PROGRAM  TO  CALCULATE  THE  2_D  BOUNOARY  LAYER  OH  A  ROTATING 

C  BOUNDARY  LAYER. 


AUS (201, XTC (4 1 


.  CER  IVX.PRINX.TE 


AUPP (25 1 


CAKMI  73  »»  Mrn  I  1*.  iai  itvrmi  I'trruji 

C0MM0NAUh5ER  CkIcuJctJchJcL,  N1.N2.NBP.DEGRF *CB . T AIJ . TH E T A . N 1 P . 

CH2P* NlMfN^HtNAf iLP^fAtPjfiCf 2PSfQl»i2f|GAM*AOELTtCIRCf|8ET # AU# AKf AH 

C, AT, AL.AKS.AKl,  AUP.AKP.AHP ,ATP  .ALP.ABP  tAUN.AKP  tAHPt ATP  t ALr . 
C^B,  3>8mJ  AUS*XZLXsT,OUXcft»UItGH,NFRST,NINTY  .NltCOP.PI 

CONSTANT  INPUT  DATA 

CATAI AUS( I). 1*1, 20  I/-. 3,- .2 5,-« 2195,-. 17647,-. Ill ,-.0526,0.0, 
C.0476* .It  .2*  .333, .5, .6, .667,. 714* .75, .8,  .857 ..889,9./ 

CATAI AKS 1 1 ) • 1*1. 20  1/1.7534* 1.7* 1.67, 1.631 *1.581 *1.543, 

Cl  51511* 

Cl ! 493. 1.474*  1.450,1.438, 1.44, 1.449, 1.455 ,1.459, 1.463 » 

Cc  AT A? Ik it  I). }»llli  1/9^000,12. COO,  7.000,  7.000,  7.000,  7.000, 

CC AT A^°fc(J  1/ 1. C/CU2/-C • 30/THETA 1/3. 93/THE TA2/C.29/CKl/3.3*8/CK 2/ 1.75 

10iCF0RMAT(//23H  INTEGRATION  STARTEO  AT  ,F9.4,4H  RAD,F10.3,8H  DEGREES! 
PI  =  DAC0S(-1 .0 ) 

CEGRF  =  180. C/PI 
ASSIGN  601  TO  NN 
5  K  *  0 
J*1 

1*0  _ 

R  E  Ac"o  AT  a'oBI  A  IN  IED-FROM-THE~  sImTl  ArI  TY^SCLUtIoNS. _ _ 

"TRiAnSWOcTBlATBETA.FnTTHATATDiLTATxlNTSG.XINTC 
I«  I«1 

IF (K.EQ.C  )  GO  TO  8 

7  IF ( G . EQ.BX I  GC  TC  10 _ 

TfImARRAYSCOUNT  THE  NUMBER  OF  DATA  ENTERED  _CERTAIN_VEL0CITY_ 

M 31= I-l”"~ _ 

” I F~a” EQUALS~Of_THE~SET~OF”oATA  OBTAInId  GOES  TO  ThI  ARRAY  SUBROUTINE 
TO  EE  FORMED  INTO  ARRAYS  _ 

IFIA.EC.0.0)  GO  TO  20 
J*J+l 
1*1 

fi  EX *P 

9  AUPPIJI  =  1.  ♦l./B 

10  K*K+ 1  _ 

INTEGRAL  VALUES  ARE  CALCULAT EO  FROM  DATA _ 

THESTR  =  XINTC  ♦  3.*8*THATA-B4B*0ELTA 
8H ( J,  II  *  B*CELTA/THATA 
ET  ( J  J  II*  -THATA*F11/B/B/B 
ELIJ.I)*  -THATA*  XINTSQ/B/B/B/B 
GK ( J , II*  THESTR/B/THATA 
E8(J,I)  =  BETA 
19  GO  TO  6 
20  L  *  20 
N*  25 

WR  IT E  ( 6, 7051  M.  K 
GO  TO  NN, (601, 602. 6031 
601  KPP  *  150 

CO  620  I  *1, LJ 
620  AUP(  I)  *  AUPP ( I » 

NIP  =  LJ 
N2P  *  KPP 
NC  *-10 
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oouu 


c 

c 

i 


c 

c 

c 


extern alFvelo!?tyl  values  wh£re  mll  velocity  IS  GREXTER  than 

5JL^0455jY*8B»  8R»BH»BTtBL,L,N,N,AKP,AHP,ATP,ALpTR*KRP7u»ABP7KC) 
■TEI6,707)  AUP 

jljjj^AMAKPnMATPlJ.n.J-i.Ul.^ltKPP) 


MR 

MR 

NR 

NR 

NR 

NR 

NR 

NR 


TE 

TE 


AUP 

RN  1 1  CIO*'”'  "  *  ' 

500  CONTINUE 

ASSIGN  602  TO  NN 
GO  TO  5 
6C2  KP  *  135 

CO  605  I  ■  l.LJ 
605  AU(I)  *  AUPPl I ) 
N1  *  LJ 


l^ililAKPU)  » (ABPCJf  n,J»lfU),I»l,KPP) 


N 


it  2  *  KP 
CC  *“10 


INVEL0ciTYALUES  m*V  FORMEO  hHERE  EXT.VEL.IS  GREATER  THAN  NALL 


RRRRTI8B»BKt8H»8T»BLf LfN»HfAK#AHtAT»AL»KtKP»LJ»ABtNC) 
GO  TO  501 
WR  IT E46>  707 )  AU 

NjjJITElfcW^HAKdMAHU.n.j.l.LJI,1*1,^! 
jjJjTIijIjCAjlAKm.IATU.n.J-i.LJI.I-i.KPI 

WR  ITE16,  ^OtI  ^AIJ  I* J“1»LJ)  *I*1»Kp) 

CONT INU|7  ^ ^AKIIIjIABNtll ,J*lyLJ)yI*l,KPI 

ASSIGN  603  TO  NN 
GO  TO  5 
603  KPN  *  75 

„  CO  610  I  ■  1.  LJ 
610  AUP I  1 1  *  AUPPl I) 

NIP  *  LJ 
N2P  *  KPM 
_ NC  =-10 

IntIgral  values  array  for’rIvIrsFflon  . 


5C1 


?nLtn*55iY,BB*  BKfBHfBTtBL,LtN,P,AKH,AHP,ATP,ALP,KfKPP,LJ,AeFfNC) 

uU  I  0  5U£ 

WRITE(6,709)  AUP 

NRITEI6I7QS)  Sum1 1 <ahmi j* 1 1  * J-i*lji .e-i.kph) 

NR  IT f ( 6*  70S J ^  AUP1 1 ^ATRI  J» 1 1 1 J*1»LJ) » 1 *1 tKPPJ 

S8  iTIS  l:?8! !  ‘  J5S"  >*  1 » .'-i.-phi 

S02  . . 

8B?B.Fg  itr  awrjypjWHffi*  ;o,nt 

30  REAUS.Tggi  UW^BEUyBETzIxCLD 


R|ACI5.7gO  UN*  BETl* BET2f  XOLO 
VRITEI6,52)  UN »BE Tl.BE T2? XOLO 
IFIUW.GT .10.0)  GO  TO  300 


LCRIT 
NEC*5 
NP AR*6 

EPCON  *  fl.OE-3 
CO  32  I  =1,5 
32  XNEWI  II  =  XOLOU) 
GO  TO  34 


CONVERGENCE  IS  ACHIEVEO  IF  CONOITICN  IS  PET. 


35  LCRIT  «  0 

IFlABSIBETl-BETNll.LT.EPCCN.ANO.ABSC  BET2-BETN2 I .IT.  EPCON) 
1LCRIT  =  1 
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THE  NEW  SEPARATION  POINTS  Of  THE  FLOW  MODEL  ARE  GIVEN  BY  THE  VALUES 
OETAINEC  FROM  THE  BOUNDARY-LAYER  CALCULATIONS. _ 

eItI=betni  _ _ _ _ _ 

"IncrI 4sl~ anglI  of  th!  lowIr  sIparatIon  point  ILC0Npnl0N.is_-i: _ 

If7aBS4BET2-bItN2T .LT.I.I-sT  GO  TO  33 
GET2*BETN2 
GO  TO  34 

33  BET2  =  BETN2  -  EPCON 

"CONDITIONS  FOR  INTEGRATION  OF  BOUNDARY  LAYER  FOR  THE  UPPER  PART 
OF  THE  CYLINCER.  _ 

~34~NFRST  *  1 

CNXT  =  l.CE-6 
CPST  -  5.0E-6 

USe”pREVIOUS  VALUES  CF  FLOW  MODEL  PARAMETERS  AS  GUESSES  FOR  THE  FLOW 
NOCEL  PARAMETERS  TO  BE  FOUND. 

~"c0*36  1=1.5 
36  XOLDI I )=XNEW( I I 

-^[[‘SUBROUTINE  TO  FINO  THE  FLOW  MODEL  PARAMETER  VALUES  GIVEN  THE 
SEPARATION  POINTS  IN  THE  PHYSICAL  PLANE  AS  INITIALLY  GUESSED  OR 
CALCULATED  FROM  THE  INTEGRAL  BOUNDARY-LAYER  EQUATIONS. _ 

CALLviLP  ARI  iln  ,  BET2,  XOLO ,  XnIw  .nIq  »NPArI 
Cl  =  XNEWt 1) 

C2  =  XNEWI2) 

AGAM  =  XNEWI3) 

AOELT  =  XNEWI 4 ) 

CIRC  =  XNEWI 5 ) 

X1C  *  0.0 
XI  =  X1D/CEGRF 
XZ  =  PI  -C.01-EPS 

IFILCRIT.EQ.il  GO  TO  95 . . . 

"cA[["sUBROUTINi~TO  CALCULATi”THi  iocl  VELOCITY  AND  ITS  DERIVATIVE 
WITH  RESPECT  TO  ANGLE ^CIVEN_A_POSITI ONONTHEC YLINDER. . 

CALL"vEL0C(ui[0Ux[xiTxZtNm . . . 

*"c7NT”is_THE"lNTiGRAL”l STARTING  FROM  X«0)  OF  THE  PART  OF  THE  LOCAL 
PRESSURE  COEFFICIENT  CORRESPONDING  TO  THE  LOCAL  HORIZONTALFORCE. _ 

CdInT~""=CsTn(X1  )-UI*UI*xI/3.G*UI*UI*Xl*Xl*Xl/10.0 _ 

CLINT  IS  THE  INTEGRAL  ISTARTING  FRCM  X*CI  OF  THE  PART  OF  THE  LOCAL 
PRESSURE  COEFFICIENT  CORRESPCNOING_TO_THE_LOCAL_VERTICAL_FORCE. _ 

cuNT*occi(xmuI*uI*xI*xI/4.c-i.o 
40  XIIND  =  0.1 

X1C  =  X1C  ♦  XIIND 
CPCOEF  =  I.0-UI*UI 

CUNT  =ClInT4X1INC*(-CPCPEFI*SIN(X1{ 

COINT  =  CCINT  ♦  X1INC*CPC0EF*CCS( XI I 

X1-X1C/0EGRF  _ 

CALL  VELCCIUI»DUX»X1»XZ»NIT I 

XZ  LAST  =  XZ 

CU-UW/UI 

CB  =  UI/IUW-UII  _ 

C  IF  THE  CYLINCER  IS  NOT  ROTATING*  A  DIFFERENT  SET  OF  INITIAL 
C  CONDITIONS  ARE  REQUIRED  FOR  THE  INTEGRAl_B-l_ECUAT I ONS. _ 

'  IMUW*UW  .LT .1 .i-fll  GOTO  50 

C  'cONCI T IOn'usIc'[n-CRdIr_TC- STip~AWAY_FRCM  THE  STAGNATION  POINT  TO 
C  A  POSITION  WHERE  THE  B-L  EQUATIONS  CAN  BE  INTEGRATED. 

- - 

IF  I CU.GT  .CU1 IGCTC  40 
CK=CK1 

THETA=THETAl 
NA  =  2 

GO  TO  55 
50  CK*I. 62574 
THETA-O. 29234 
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5;3puI*?o  kutner 

55  CONTINUE 


52Sc.I.II9Ss.ilSEP,Sy5I!l!§  integration  of  the  boundary  layer  for  both 

THE  UPPER  ANCJ.CHER  PART  OF  THE  CYUNOER. 

fco’xTII~*xI""- 

X 2)-THETA«TFETA 
X(2)»CK*CK*xm 
X(A)  -  CLINT 
X ( 5 )  >  CCINT 
70  CM  AX  «  2.0E-3 
ER-  5.0E-6 
NINTY  *  1 
PS«1.0 
HR ITE( 6. 150 1 

CALL  SUBROUTINE  TO  INTEGRATE  BOUNDARY-LAYER  EQUATIONS 

CALL  KUTHErTcNXTi CPST«OMAXf  5t  X » OXtDERI VX «ER  •  1 1 H t0 «PS «XPVf PR  I  NX « 
CXTCt  A  «NTS. TERMX) 

GO  10(75,76.17.781, NTS 

75  UR  IT E ( 6 •  1 10 1  NITCOM 

IF (NFRST .EQ.l )  GO  TO  73 
EETN2  «  X( 1 1-0.C05 
GO  TO  35 

73  8ETN1  *  X(l»  4C.01 
GO  TO  85 

77  UR  ITE(6. Ill) 

GO  TO  35 

78  UR ITE( 6, 112 ) 

GO  TO  30 

76  UR ITE( 6.123 )  X( 1  ),X( A) ,X( 5) 

IF (NFRSi «EQ • 1 ) GO  TO  85 
IFINFRST.EQ.2)  GETN2  »  X(l) 

GO  TO  35 

COMPU TE~ INPUT~FOR~LCuiR-PAR T~CF~CYL INDlRr__siI~COPplNTS~FOR_ I NPUT 
T 0  UPPER  PART  OF  THE  CYLIJOER. 

85~NFRST*2" 

EETM  =  XIII 
CLUPPR  =  X ( A ) 

COUPPR  =  X ( 5  > 

X1C  a  41.0E-5 

XI  =  X10/0EGRF 

XZ  =  -PUG. 01  4EPS 

CALL  VEL0C(UI,CUX.X1>XZ,NIT ) 

CLINT  *OCOS(Xl)4UlnJl*xi*Xl/A.C-l.C 
CCINT  »CSIN(X1)-UI4UI*X1/3.04UI*UI*X1*X1*X1/10.0 
88  X1C  »  X1C-X1INC 
CPCOEF  *  1.0-UI*UI 
XI  INC  =  -XI IND/CEGRF 
CLINT  =CL  INT  4X1 INC*( -CPCOEF )♦ SIN ( XI ) 

COINT  »  CCINT  4  X1INC*CPCCEF*CCS(X1) 

X1«X1C/0EGRF 

CALL  VELOCI U ItOUX.Xlf XZf NIT ) 

XZ  LAST  a  XZ 

CU*UU/UI 

CB  a  UI/(UU-UI) 

CNXT  *  -1.0E-6 

CPST  =  -5.0E-6 

IF (UH*UU .LT.1.CE-6I  GO  TO  50 

IF(CU.LT.CU2.0R.CU.GT.l.E-7)  GO  TO  88 

CK*CK2 

THETAaTHETA2 

NA»1 

GO  TO  55 

T  Hl"cONvlRGlNci"cR  IT  Eiu0N~F0iTTHl”  sIPARATION~  P0INT5’I  S~S  AT  IS  Flic” 

THE  COEFFICIENTS  OF  LIFT  AND  DRAG  ARE  CALCULATED  ANC  URITTEN  OUT. 

sI  coIfl  »  zIo4PI*circ*cos( alphT* 

HR  ITE(6*  710)  EPCON 

CO  EFLF=( CLUPPR-XI A )-( 1.0-UIC0M*LIC0M) * (COS ( BE TNI ) -COS (BETN2 ) 

C))/  2.0 

CO  EF0Fa( CCUPPR-XI 5  )♦( 1 .0-UICOMUICGM)*(  SIN  (BETN2 1-SIN  (BETN1 ) ) 1/2.0 
UR ITE(6* 715)  COEFL 
UR ITE( 6. 720 )  COEFLF 
UR ITEf 6»  725)  CCEFDF 
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GO  TO  30 

52  F0RMTUH1.30*  SPIN  RATE  AND  INITIAL  GUESSES  /11X,2HUW.6X, 

C11HSEP.  ANGLE  1, 5X, 1QHSEP .ANGLE2, 8X.7HSGURCE1 , 6X .7HS0URCE2 , 

C5X . 10HSORC  ANGL1, 5X, IOHSORC  ANGL2.$X,12HCIRC.  v£lCC./8F15.4 ) 
i30  FORMATION  ERROR  RETURN. ..ITERATIONS  IN  VELOC  SUBROUTINE  HAS  EXCEE 
CCEC  ,  110) 

111  F0RMATI23H  K  EXCEEDED  UPPER  UNIT) 

112  FORMATI 23H  K  EXCEEDED  LOWER  UNIT) 

123  FORMATI 13F  SEPARATED  AT  ,F12.4,8H  RADIANS  ,/31H  LIFT  COEFFICIENT  C 
CONTRIBUTION.  ,F12.5,/31H  DRAG  COEFFICIENT  CONTRIBUTION.  , 

CF1/.5) 

150  FORMATI 1F0, 6X, 1HX.8X.5HTHETA, 7X,1HK»6X»5HUW/UE ,5X,5HDELTA,6X» 
700CFORiAf(8flo  *6X. 5H0U/DX.8X.2HCL.8X.2HC0.8X.2HCP/ ) 

701  F0PMATI1F1,(7X»16F7.3)) 

7C2  FORMATI 1F0, I 17F 7.3) ) 

7C3  FOPMATI1HO.I 15F8.4 1 1 
7C4  FORMATI 1F0, I F8.4, 17F7.3) ) 

705  FORMATIIHG,!  2  C 1 1 0  I ) 

706  F0RMAT(1F1,(6X,14F8.4) ) 

7C7  FORMATI  1U»  I  8X.17F7.3)) 

708  FORMAT! lHO. I 15F8.4) ) 

7C9  F0RMATI1H1.I8X.15F8.4) ) 

710  F0RMATI//26H  ANGLE  INCREMENT  LESS  THAN  ,  E10.4,1CH  SATISFIED  ) 

712  FORMAT  ( 24F  SEPARATION  VALUE  OF  H...F10.5  ) 

715  FORMATI/21H  LIFT  COEFFICIENT  IS.  .E12.5) 

720  FORMATI/22H  LIFT  COEFFICIENT  BY  INTEGRATION/  E12.5) 

725  F0RMATI/22H  CRAG  COEFFICIENT  BY  INTEGRATION/  E12.5) 

3C0  STOP 
ENC 

C  * 4 **** ****** * ******************************* *************  ********** 

C  SUBROUTINE  TO  CALCULATE  THE  DERIVATIVES  FOR  DIFF.  ECUATIONS. 


SUEROUTINE  OERIVXtX.DX) 

CIMENSION  XI5J.OXI5) 

C I  MENS  ION  AUI17),AKS!20).AKH2C).AUS!20),  AK 1 1 35 )  ,AH 1 17 .135 )  , 

XZ  =  XZL AST 
X1=X( 1) 

CALL  VELOCI U I.CUX.Xl.XZ.NIT) 

CU  =  UW/UI 
Cfi  =  UI/IUW-UI) 

CKSG  =  XI 3 )/X( 2  ) 

CKSCA  *  ABStCKSQ) 

CK  =  SQRTICKSQA) 

SELECT  SET  OF  ARRAYS  FOR  PARTICULAR  VELOCITY  RATIO  ANC  OBTAIN 
H.TfL.BcTA  FCR  VALUE  OF  K  KNOWN. 

GO  TO  (1,2,3), NA  ”  . 


:  ch 

=  TDINTICU 

»CK 

V 

AUM 

»N1M 

9 

AKM 

9 

N2M 

9 

AHM 

9 

N1M 

f 

*■> 

c 

9 

2) 

CT 

=  TCINTI CU 

,CK 

9 

AUM 

,N1M 

9 

AKM 

9 

N2M 

9 

ATM 

9 

N1M 

9 

2 

9 

2) 

CL 

GO 

=  TCINTICU 
TO  4 

,  CK 

t 

AUM 

,N1M 

9 

AKM 

9 

N2M 

9 

ALM 

f 

MM 

9 

2 

f 

2) 

2  CH 

=  TCINTICU 

,  CK 

V 

AUP 

,  NIP 

9 

AKP 

9 

N2P 

V 

AHP 

9 

MP 

V 

2 

9 

2) 

CT 

=  TCINTICU 

,CK 

9 

AUP 

,N1P 

9 

AKP 

9 

N2P 

9 

ATP 

V 

NIP 

9 

2 

9 

2) 

CL 

GO 

=  TCINTICU 
TO  4 

,CK 

9 

AUP 

,N1P 

9 

AKP 

9 

N2P 

» 

ALP 

9 

MP 

9 

2 

9 

2) 

3  CM  «  TCINTI 
CL  =  TCINTI 
CT  =  TCINTI 
GO  TO  4 

4  CX(l)  =  1.0 


CU,CK, AU,N1,AK,N2,AH,N1,2,2) 
CU,CK,AU,N1,AK,N2,AL,N1,2,2) 
CU,CK,AU,N1,AK,N2,AT,N1,2,2) 


INTEGRAL  MOMENTUM  EQUATION 
CXI2)  =  I4.*CT-0UX*XI2)*(CH42.C)*2.0)/UI 


INTEGRAL  ENERGY  EQUATION 


CXI3)  =  I8.*CK*(CL*CU*CT)-6.*X(3)*DUX)/UI 
CPCOEF  =  l.U-  UI*UI 
CXI4I  =  -CPC0EF*SIN(X1) 

CX  c  5  )  =  CPCOEF*CCS( XI I 
X2  A  =  ABSIXI 2 1 ) 

THETA  =  SORT  I X2A )*ABS(CL ) /CL 
TAU  =  U I ♦CT/THET  A 
CUXCN  =  CUX 
UICCM  =  UI 
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GO  TO  (5t 6*7 ItNINTY 

5  XII AST  ®  XZ 
RETURN 

6  XZ  LAST  *  XZ-C.02 
RETURN 

7  XZLAST  *  XZ+  0.01 
RETURN 

C  ****SUBROUTINE  TERMXI X, 0X,XTC,XP\n _ 

£  "SUBROUTINE* FOR^iRMINATING  ThI  InTEGRaIcN _ 

HSInI  Jon  5&(  AK11 201  »AjJS( 20,  »  AK(  135,  •  AH (17^135 )j 

CAKM(  75 )«  At*H(  14» 75 ,*ATM( 14»75, *ALM( 14»75,  .ALM14, 

C0PP0MUSE  HA  IN  I  _ 

VALui~OF~NA~is”siT~TG”siLECT  RIGHT  ARRAYS_OF_INTEGR AL_VALUES. _ 

na“T 

IFICB.GT .0.0 )  NA  =  2 

IFICB.LE.-. 959999)  NA  =_3 _ 

obtaik  separaticn  valueicksT  of  k  PCR^ociTvjATio.-^: _ _ _ 

. C  AL  i”cv  C INT  ( CUTCK  sT  AUS ,  AKS.NBP.il . . 

IF  ANY  ELEMENT  OF  XTC  CHANGES  SIGN*  INTEGRATION^^  STOPPED. _ 

XT  cl  u”®”  29  •  C*"“”fL  OAtTnI  TCOM ) 

1  XTC(i!S=KBET22  -XII,  .0.001 
EET22  =  EPS-PI  ♦  2.*ALPH 
XTCf 4 »  =  1.0  , 

XPV  =  XIII  *  OEGRF 
RETURN 

r  *****!»******>»***♦****♦**♦**’*********♦♦♦♦*****♦♦*♦******************** 

SUBROUTINE  PRINXIX.CX] _ _ 

PRINTING  SUBROUT INE  »♦♦*«**♦***♦*♦*»*♦*♦*♦♦***♦♦♦♦*♦♦♦♦♦***♦ 

1141 

COPPONCUSE  MAIN  I 
FBCOM  *  CHLAST 
XBCOM  *  XCLAST 
X1C=X(1,*CEGRF  . 

cpcoef  *-cxm/siN(xm, 

WR  IT E, 6»  10* THETA tCK.CUt DELTA »TAU*CL»UI COP  .OUXCNi  XI4),  X(5,  , 

CIF(ABS(XlC).GE.e9.C.AN0.NFRST.|Q.l{  NINTY  =2 
IF(ABS(X1C,.GE.89.0.ANO.NFRST.EQ.2I  NINTY  *3 
XOLAST  =  XC1) 

CHLAST  =  CH 
1  FORMAT! 12F10 .5  ) 

RETURN 

r  MM*Mt*«*mM***MM******«M«*M**M**«MM*MOM»MMMM<M 
SUBROUTINE  VELOCI U I,DUX,X,XZ_ .NIT] _ 

C  _SUBR0UT IN E~ To” COMPUTE* ThI” IN vl SCIo” vIlOcTtY  AT  A  POINT  ON  THE  CYLINC* 
C  ER  IN  THE  PHYSICAL  PLANE.  _ 

COMMONIUSE  MAIN) 

NI  T=0 

CALP’»*CCCS( ALPH) 

. 

c  - - 


onoo  oooooo  o  ooo  ooo  oooo 


C  ITERATIVE  METHOC  OF  NEWTON  TC  DETERMINE  CORRESPONDING  ANGLE  IN 
C  UNTRANSFORMEC  PLANE. 

C - - - - - - 

1  XZC=XZ 

COSX  =  CCCSI XZO  +  EPS I 

ThI”c0RRlcT~RlLATl0N"ls’G7viN**BiTWEEN  CORRESPONDING  POINTS  IN  THE 
TWO  PLANES  WHEN  FX  IS  ZERO. 

Fx”=”7I»0-CALP*COSx7*DsTN7xZo”+Ipi7/7 Al-CCSxI~SINX 
FPX=(Al*(COSX-CALP*<2.0*CCSX*CCSX-1.0H-l.C*CALP*CQSXM3)/ 

C(  A1-C0SX )<*2 
XZ*XZO-FX/FPX 
NIT=N  IT*1 

IF  IN  IT .EC .30 )  GO  TO  2 
IF  I ABSIXZ-XZG ) »GT .PR EC IGOTC  1 

2  SINXZ  =  CSIN(XZ*EPS) 

COSX  =CC0S(XZ4EPS) 

GF=2.0*SIN(XZ)-Q1*CCTAN(0.5*(XZ-AGAM))-C2*CSTAN(0.5*(XZ-ADELT)) 
C*CIRC 

GS  =  1 .  0-2.C*CALP*CCSX*CALP*CALP 
G=0F*GS 

H=2.0MCALP-CCSX) 

_  —  —  mm  —  m*  —  —  m.  —  m,  —  ——  —  ——  —  —  —  —  mm  mm~***  —  m*m*  —  mm  —  —  mm  — —  m*  ————•*  -m  —  —  mt  m,  —  ———————————————  —  —  f* 

THE  ECGE  VELOCITY  FOR  THE  GIVEN  PCINT  CN  THE  CYLINDER. 

UI  =  G/H 

GP  =  (  2  .♦OCOSI  XZ)-».5*Q1/(DSIN(.  5* (XZ-AGAM)  l**2)«-.5*C2/(DSIM.5* 
C(XZ-ACELTn**2n*GS*GF*(2.0*CALP*SINXZ) 

HP  =  2 . 0  *  S INXZ 

XNUME  =  !  Al-COSX)*! COSX+CALP* II  SINXZ**2)- (COS XM2||| -SINXZ* 

CIS  INXZ-CALP*COSX*S INXZ I 
XDENO  =  ( Al-CCSX )**2 
YD  EN  =  XNUME/XCENC 


THE  SPATIAL  CERIVATIVE  CF  THE  VELOCITY  AT  THE 


THE  6-L. 


CUX- ( H*GP-G*FP )*OCGS (X-EPS ) /H/H/YDEN 

MTCOM  a  NIT 

RETURN 

ENC 

4*^***^ ************************************************************** 

SUBROUTINE  VELPAR(BET1,BET2,XCL0,XNEW,NEQ,NPAR) 

SUeROUT INt_TC~CALCUL ATl~THE  PARAMETERS  FOR  THE  FLOW  MODEL  GIVEN  THE 
SEPARATION  POINTS  IN  THE  PHYSICAL  PLANE.  THIS  SUBROUTINE  USES 
CTHER  SUBROUTINES  TO  OBTAIN  THESE  PARAMETERS  BY  A  NEWTCN-R APHESON 
SCHEME  FOR  SOLVING  A  SYSTEM  CF  NONLINEAR  EQUATIONS. 

c7mENSION~a7sj6 )»p7i#6 1 tQoTs.tT *TrT 5 1  * JC ( 5 ) *cTi *6 > ,X(5) 

CIMENSION  XOLCINEG ). XNEWINEQl  ,  . 

C I  MENS I^N  AU(17).AKS(20).AK1(2C),ALS(20),  AK 1 135 ) .AH  1 17 .135 ) , 
CAKP(160),AHP(17,l6C»tATPI17,lfcC!.AlP!17,16C!  ,AUP(17) ,ABPI 17,160) , 
CAT  1 17, 135  I. AUI 17,  135 ), ABI 17, 135 1, ABM (14, 75 > , 
CAKM(75I.AHM(lAt75).ATM(lAf75)fALM(lA,75f ,ALM(1<*> 

CIMENSION  XNEWBI 5  I 
COMMONIUSE  MAIN) 

EXTERNAL  JACOB 
FI  =  CACOS(-l.O) 

EPS=  I  BET  1  *  BET2 ) *0.5 
NM AX  =  NEQ 
N  =  NEQ 
M  =  NPAR 
L=  0 

UMAX  =  1000 
EPS!  =  1.0  E—  7 
EPS!  =  l.CE-5 
S I  =  i  .0 
CS  =  ! .  0 
SF=5.0 
EPS2= 1 .OE-2 
CO  10  1=1, N 
CO  5  J=l, M 
PI  I,J)=0.0 
00(1, J)=0.G 
5  C(  I, J  )  =  0.0 
10  CONTINUE 

ALPHISHALFTHE  ANGLE  BETWEEN  THE  TWC  SEPARATION  POINTS  IN  THE 
UNTRANSFORMEC  PLANE. 

~’aLPF=C.5*(p7-BET1«EPS) 


nnoo  n  o  onn 


YN  «  0.1 
20  YN-0.5*YN 
YNN  *  VN 

25  AINC  -  PI*YNN/ieO.C 
ALPINC  »  ALPH  -AINC 


CALL  SUBROUTINE  THAT  SOLVES  A  SYSTEM  OF  NONLINEAR  EQUATIONS. 


ttKAX,|fljJC,NyM,L,ITNAX,SI,OS,SF,EPSl, 


*■108)  ITMAX 

_ it  105)  EPS2 

IFUERRI.LE.C.CR.IERR2.LE.0)  IJK  >! 

WRITE! 6, 107)  !(XN£V(I),I«1,N),YNN) 

IFIYN.LT. 6. OE-2.ANC.YNN.LT.C.0)  60  TO  50 
CO  30  1-1,5 
30  XNEWB! I ) -XNEW! 1 ) 

IF(YNN.LT.O.C)  GO  TO  20 
YNN  =  -YNN 
GO  TO  25 
50  CO  60  1-1,5 

60  XNEW!  I  )-0. 5* (XNEW! I) ♦XNEWBI  I) ) 

YNN  -  0.0 

COEFL  *  2.0*PI*XNEW(5)*C0S(4LPH) 

WR  TE(6,1C$)  BET 1,BET2* ALPH, EPS  , COEFL 
WR  T  E ( 6 , 110 ) 

WR  TE(6,107)((XNEW(I),I-1,N),YNN) 

2C0  RETURN 

WRITE!6.106)((Q( 1,J)VJ»1,M),I»1,N) 

102  F0RMAT(8F  I  ERR 1  =,I4//) 

103  F0RPATI62F  ERROR  RETURN - JACCBIAN  MATRIX  IS  SINGULAR  OR  ILLCONCI 

CTIONEC  ///) 

FORMAT! 8^  JERR2  •=  tI4//) 

FORMAT! 38H  ERROR  RETURN. . .STEP  SIZE  IS  LESS  THAN  ,016.8///) 

106  FORMAT! //20H  FINAL  Q  PARAMETERS.  //(5D16.8)) 

107  F0RMATI//29H  FINAL  SOLUTION  FCR  THE  ROOT. /8X,2HQ1 ,14X,2HC2.12X, 
C5HGAMMA, 1 IX.  5HCELTA, BX. 1 1HC IRCULATICN, 5X.11H ANGLE  lNCR./<6Fl6.8l ) 

3  F0RMATI//58H  ERROR  RETURN .. .MAXIMUM  NUMBER  OF  ITERATIONS  HAS  EXCEI 


104 

105 


1C8 


CCOTAN! X ) 
COTIX ) 


CACCSIX) 

ACOS(X) 


109  F0RMAtl/56H  PHYSICAL  PLANE  ANO  PRIMARY  TRANSFORMED  PLANE  PARAMETER 
CS  »//8X,6H  BETA1, 10X,5HBETA2, 10X,5HALPHA ,9X ,7HEPS IL0N,7X,9HLIFT  CO 
1EF  ./5F15.5//  ) 

110  F0RMATI//24H  BY  LINEAR  INTERPCLATI ON) 

RETURN 

ENC 

** ,44* ***************************************** ******4***** ******* 4* 4 

FUNCTION  COTAN! X) 

COTAN  -  COT! X ) 

RETURN 
ENC 

FUNCTION 
CCOTAN 
RETURN 
ENC 

FUNCTION 
CACCS  » 

RETURN 
ENC 

*** ****************************4************************************** 
SUBROUTINE  JACCB(Q,NMAX, aNEW,N,M,A) 

SUBROUTINE  TO  COMPUTE  THE  JACOBIAN  OF  A  NONLINEAR  SYSTEM  OF  ECUATIONS 
CIVEN  THE  TRIAL  PARAMETERS  ASSOCIATED  WITH  THE  FLOW  MODEL. __ _____ 

CImInS ION  AU ! m , AKS ! 2cT . AK1 1 2C # ,  At S ( 20 ) ,  AK ( 1 35 ) . AH (  17 .135 ) , 

CAKP! 160 ) , AHPi 17, 160), A TP (17, 160) , ALP! 17, I 60 ) , AUP (171 ,ABP( 17,160) , 
CAT!17,135)*AL!17»135),AB(17,135)tABM(14«75)» 

CAKM! 75), AHM! 14. 75 ) , ATM! 14, 75) ,ALM( 14.75 J , AUM (14) 

CIMENSION  Q(NMAX,MJ,XNEW(N) ,A(NMAX,M) 

COMMONJUSE  MAIN) 

CALP  -OCOS! ALPH ) 

CALPIN  -CCOS! ALP  INC ) 

-  -CS IN! ALP  INC  ) 

-CCOTAN! 0. j*(ALPH-EPS-XNEW(3)>) 

=  CAMEM3 

=  CC0TAN(0.5*(ALPH-EPS-XNEW(4) ) ) 

-  CAMEM4 

-CSIN(0.5*(ALPH-EPS-XNEW(3)>) 

-  0.5*XNCW(  D/SAMEM3/SAMEM3 
=CSIN(0.5*(ALPH-EPS-XNEW(4) 1 1 

-  0 .5*XNEW( 2 1/SAMEM4/SAMEM4 
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SALPIN 
CAMEM3 
A!  1,1) 
CAMEM4 
A!  1,2) 
SAMEM3 
A!  1,3) 
SAMEM4 
A!  1,4) 


nno  noon  onn  ooo 


A ( 1,51  =  -l.C  _ 

A?  1, 6  l=-lxNE(i<  1  l*CAMEH3*XNEht  2  l*DAMEMA-2.G*SAME-XKEk(5 )  ) 

CApIp  3  »CC0TAN(0.5*(ALPH*EPS*XNEk(3m 

CAPlpA  *CC0TAN(C.5*<ALPH*EPS*XNE»((Am 

SAPEP3  *DSIN( 0» 5*<  ALPH*EPS*XNEh43n  ) 

A ( 2,3)  *  -0.5*XNEkil  I/SAPEP3/SAPEP3 
SAPEPA  =  CS  IN( 0.5*( ALPH+EPS+XNEV ( Al ) ) 

A(  2, A )  *  -0.5AXNEK2I/SAPEPA/SAPEPA 

A? 2io l=-( xSeS^I j*CAPEP3*XNEtal 2  )*CAPEPA-2. 0*SAPE+  XN'EM  5 ) I 

ACM&M3  =  CC0TAN(C.5*(ALPINC-EPS-XNE|.  3 

ACPEP3  =CCOTAN(C.5*ULPlNOEPS*XNEh<3m 

A(  2.1J  =  -ACHEM3*ACPEP3 

ACPEPA  =CC0TAN( 0 • 5 * < ALPINC+EP  S* XNE K I A I ) ) 

ACMEMA  =CCCT  An! ct 5*( ALPINC-EP  S-XNEb( A I ) ) 

A ( 2.2  *  =  ACPEPA-ACMEMA  , 

SCP^PS  =CSIN(0.5*(ALPINC*EPS+XNEW  3 

SCPEM3  =  CS  IN  ( C.5*(  ALPINC-EP  S-XNEW(  3  Ml 

AI2.il  =-0.5AXNEUI  11*1 1.0/SCNEM3**2+1.0/SCPEP3**2I 

SCP^PA  =CS  INIC .5*(  ALPINC*EPS+XNEVI(A I J  j 

SCHEMA  =CSIN(0.5*(  ALP  INC- EPS- XNEtal A  III 

A ( 3, A)  =  -0 • 5*XNEW(  2  )*( 1.0/SCNEMA**2+1.C/SCPEPA**2) 

At  2,51  =  2.0 

SACPE  =CSIN(ALPINC*£PSj 

A?  3ii6 1=-?2 .0*1  SACME-IaCPE  I-XNEM1)*(ACMEM3-ACPEP3)-XNEW(2)*(ACMEMA 

CCCMEM3A=CCOS(C°5*l  ALPINC-EPS-XNEW(3I) ) 

C  =  0:5*(1.0-2.C*CALP*CALPIN4CALP**2> 

C  =  SALPIN*(CALP**2-1.0 ) /(CALP-CALPIN I 
A ( A, 1 )  =  C/SCMEM3**2  -  0*ACMEM3 

A(a!3)  =  C*XNEHI1)4CCMEM3/SCMEM3**3  -D*XNEW (1) /SCMEM3**2/2.0 

SK5«  :cES§Aii?j{^A?5^s?S15i5,,^)cNE.c2,/scvE^-2/2.o 

At  A* 5 1  =  «0 

A?At!6)*-2^0*C*I2l!lc4CACME4C.5*(XNEWm/SCMEM3**24XNEhC2J/SCPEMAM2 

C)l-fc*I2.3*SACME-tXNE«(ll*ACMEM34XNEM(2>*ACMEPA-XNEW(5m 

nmwstizmsi  ,tll 

Sf \  §i  :cS2USfs*{i^iJ!?lcil^l2yj.5.D.x«.  a  i  /scpep3..2 
SfjiJi  :c'8ijiSf5?i?iEi^I?ic*pTOJ{.5.i).*Kew(2./ScpEp3..2 

«fE?M~(2*o'c«l2C.C*C»CPEx0.5P(XNEEm/SCPEPJ.J2*XPm2[/SCPEP4*>2 
Cl  |  -C4(-2.0*SACPE*XNEW(1)*ACPEP3  4XNEW ( 2 ) *ACPEP4  ♦  XNEW(5))I 
RETURN 
ENC 

*♦*♦**************♦*******************♦*♦****♦♦♦****’*♦*****=•***♦***** 
SUeROUTINE  ARRAY1BB* CK,CHfCTtCL,LfN,M,AK tAH^AT, AL,I< ,KP.LJ»AE,NC)_ 

SUBROUTINE  TO  FCRM  KtH,T, L.BETA,  ARRAYS. _ _ _ __ 

C  I  MENS  ION  CKtL»N)  »CHtL  ijj  j  tCT(L,N)  ,  AK  (KP)  ,CLIL » N  !) , 

CAH  (LJ.KP ) , ATILJ.KP )  , ALILJ.KP  )  ,MtU  ,BB ( L • N ) , AB ( L J ,KP J 
0  IHERs  ioft  AKXtko).  ADt3c6),CKJ(25I,CK.Jl(25j,CKJ2(25), 

C  CEJ(25I,CTJ(25),CLJI25)  ,CK J3 ( 25> ,CB J (25 1 ,ACS ( 25 J 

11  =  0  _ _ _ _ _ _ _ _ 

PLACE  ALL  TEE  K  VALUES  INTO  A  i-C  ARRAY. _ 

CO  10  J=i,LJ 

MJ  =  M{  J  I  _ _ _ 

~PLACE’ALl"K's’wiTH’sAMi”vELOcTTY”RATlc  INTO  1-0  ARRAY  AND  SORT 
IN  INCREASING  ORCER.  _ _ _ _ 

coTI’i.mj 
3  C K  J (  I )  =  CKU.II 

CALL  SGRTXYICKJ.ADS.MJ^ . . . . . . 

CONST PUCT  ARRAY  OF  cTEFiRENcis  BETREiN  SUCCEEOI NG^K  VALUES. . 

ACS(l)  =  c.o 
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4  AOS?  I)  3'CKjfl)  -  CKJ(I-l) 


PLAci”sMALLER  ARRAYS  INTO  LARGER  1-0  ARRAYS  OF  H  VALUES. 


CO  5  I»1 
II  -  II 
A0( II )  3 
S  AKX( III 
10  CONTINUE 


5  1*1* MJ 
-  II  ♦  I 
II)  *  ADS  I 
(II)  -  CKJ 


Sill 
KJ(  I) 


SURT  LARGE  ARRAYS  SO  THAT_K_IS_IN_INCREASING_OROER _ 

CALL~SORTXY( AKX.AO.tO _ _ 

" SSW iT iSr5a5i!lH?5 TSS5ili I? SEVlH!HESCtOsIoElSTlftv»;  ||¥«l"i,'G 

SUCCE EC  ING  K^G^Ie’tHE  VELcIlt"  ffllO  U  HELD  COKSTMIT. . 


J  3  2 

AK II)  =  AKX(l) 

IF  I  AC (  II.LT.l.E-6)  GO  TO  15 
AKXM  «  AKX(I)  -  ADI  I ) 

IF( AKXH.EQ.AKf J-l  j  GO  TO  15 
IF(AKXM.GT.AK(J-1))  GO  TO  13 
GO  TO  20 
AK  ( J )  3  AKX(I-l) 

GO  TO  19 

AK ( J )  *  AKX(I) 

J  *  J  ♦  1 


20  CONTINUE 


"”InLARGE**NO  .”of”nON-ZERO  ELEMENTS  TO  INCLUDE  SOME  LARGER 
VALUBS  OF  K. _ 

xT^Cu”"lAKU-n-AK(  J-2n/lC0. 

R I  =  1.0 
JP1  *  J«1 

AK(I)  *  AK(J-1)*RI*XINCU 
RI»RI  ♦  1.0 

29  CONTINUE  _ 

FORM  2-C  ARRAYS.  _ _ _ 

36  CO  70  J  *1* L J 
MJ»M< J) 

37  CO  40  I3 1*MJ 
CK J(  I  )  =  CKIJ.I 
CKJKII  3  CKJ  J 
CKJ2(  I)  3  CKJ | 

CKJ2U)  3  CKJ1 
CH J( I j  3  CH  J. 

cGjii  :  : 

*°  CALL  SORTXYICKJ.CHJ.MJ) 

Bit  BHK  Si  :gj:H 

CALL  S0RTXY(CKJ3tCBJ»MJ) 

ASSIGN  46  TO  KGT 

45  CO  60  I  3  l.KP 

GO  TO  KGTi (46*60) 

46  IP 6  3  I  ♦  6 

TO  6C 

51  }f(AK(I-6).GE.CKJ(MJ))  ASSIGN  6C  TO  KGT 

1  folfflBiftlgHAUgU:  K8:Si:8:!feI! 

AHIJ.I)  3  AH  I 
AT ( J  *  II  3  AT 
AL ( J*  1 1  3  All 
AB( J* I )  3  ABI 
60  CONTINUE 
70  CONTINUE 
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oonnn 


101  FORMAT 4 1MJ,  ( 
ICO  FORMAT! 8F10.5 

102  F0RMAT(lMj,( j 


105 

tC0 


FORMAT! 1F0,( 10F12.5) ) 
‘4F5,3|) 


P0RMATI1H 
RETURN 
ENC 


0110) 


SUBROUTINE  0ISC0T  (XA, ZA,  TABX, TABY.TABZ ,NC ,NY , NZ , ANS ) 

subroutTne  for  lagrangian  interpolation  and  extrapolation! 


08-01-68  SUBROUTINE  REVISEO  08-01-68  **#*44*** 
THIS. SUBROUTINE  ARE  ONLY  DUMMY  DIMENSIONS. 


ic 


15 

20 

25 


25,25,20 


35 

40 


45 

50 
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***  DOCUMENT  DATE 
THE  DIMENSIONS  IN 
CIMENSIQN  TA8X 
DIMENSION  TABX 
CALL  UNS  (NC. 

IF  ( NZ-1 )  5,5,10 

CALL  OISSER  < XA.TABX4 1 ), 1,NY, IOX.NN) 

NNNS I CX* 1 

CALL  LAGRAN  ( XA, TABX (NN) , TA8YINN ) ,NNN, ANS) 

ZARG=ZA 
IP  IX  =  ICX  ♦  1 
IP 1Z  =  ICZ  ♦  1 
IF  ( IA)  15.25,15 
IF  4 Z ARG-TABZ f NZ ) > 

ZAPG  =  TABZ(NZ  I 
CALL  CISSER  (ZARG.TAPZU),  1,NZ, IDZ.NPZ) 

NX  *NY/NZ 
NPZL=NPZ«IOZ 
1=1 

IF  (  IMS)  3C.20.40 

NPY( I)=< JJ-1)*NX+NPX<1) 

NPX ( I )=NPX( 1 ) 

1=  HI 
GOTO  50 

CO  45  JJ=NPZ,NPZL 
IS=(JJ-l)*NX4l 

CfSSER  4XA»TA8XI l».IS.MX»IOX.NP*IIM 

Nr  Y 111  -Nr X I  I ) 

Is  HI 

CO  55  LL  =  1,  IP1Z 
NLOC=NPX ( LL ) 

NLOCY=NPY(LL ) 

CALL  L AGR AN( XA.TABXINLOC ),TABY(NL0CY),IP1X,YY(LL)) 


70  return*gr<in  ,z*rg*tabz<np*,»VV(1,»ip1hans? 

ENC 

*4 4*** *************************************************** ************ 
SUEROUTINE  UNS  ( IC, IA, IOX, IDZ , IMS ) 

IF  (  IC  )  5,  5, 10 

5  IMS=1 
NC=-  IC 
GOTO  15 
10  IM  S  =  0 
NC  =  IC 

15  IF  ( NC-10C  )  20,25,25 

20  I A  =  0 
GOTO  33 
25  I A  =  1 

NC=NC-10C 
20  IDX=NC/1Q 
I0Z=NC-ICX*1C 
RETURN 
ENC 

4 4*****4*****#*#************4**********************4***#*****# ******* 
SUEROUTINE  CISSER  ( XA, TAB, I ,NX, ID, NPX ) 


CISC3T 


TAE(2j 


CIMENSION  _ 

DIMENSION  TAE( 2  I 
NPT=IC*1 
NPE=NPT/2 
NPU=NPT-NPB 
IF  (NX-NPT)  10,5,10 
5  NPX= I 
RETURN 

10  NLOW=  HNPE 

NUPP=HNX-(NPU*1) 

CO  15  1 1  =NLOh,NUPP 
NLCC=  1 1 

IF  C  TAB!  II)-XA)  15.20.2C 
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CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

CISCOT 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

UNS 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

CISSER 


oonr»n 


15  CONTINUE 

NPX«NUPP-NPB4l 

RETURN 

20  NL-NL0C-NP8 
NU-NL+IO 
CO  25  JJ’NL'NU 
NO  IS* JJ 

IP  (T  A8( JJ )-TA8( JJ+1 I I  25.3C.25 

25  CONTINUE 
NPX*Nl 
RETURN 

20  IF  (TAB(NCIS)-XA)  AO. 35.35 
25  NPX*NCIS-IO 
RETURN 

AO  NPX=NCIS*1 
RETURN 
ENC 


SUBROUTINE  LAGRAN  (XA.X.Y.N.ANS I 
CIHENSION  XC  2 1 .V( 2 1 
SUMO.O 
CO  3  1=1, N 

PROC=Y(It 


J*1»N 

1-xi J» 


A«X( 1 1— X f  J >  , 

IP  (At  1.2.1 
1  Ek (XA-XC  J I  I/A 


nm 

£  Lf* 


CISSER 

CISSER 

CISSER 

CISSER 

CISSER 

LAGRAN 

LAGRAN 

LAGRAN 


LAGRAN 

LAGRAN 


PRCC=PR0C*6 


2  CONTINUE  „ 

3  SUP=SUH4PR0D 
ANS=SUH 
RETURN 

iEKSa*tlt* . . 


CIMiNsioN'plNMAX.NT.QOiNMAX.MT.ATNPAX.Ht .IrInI  ,JC(M ,C(NMX,F) .XOLNSIMEQ 
1C(N).XNEW(N)»X(Ni  _ 


LAGRAN 

LAGRAN 

LAGRAN 

LAGRAN 

LAGRAN 


NS [MEG 
KS I M  EC 
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noo  oooo  oooooo  ooooo  oooo  oooo  oooo  ooooooo  .*~'Ooo  ooooo 


PROGRAM  IN  IT  III IZAT ION. 


r . - _ _ NS  IHEQ 

I E  PR  1=0  KCIWFC 

I E  PR2*Q  t5|8|S 

IF  (L.EU.l)  S * S I  NSINEC 

_ _ ns  in  eg 

TRANSFER  THE  OLC  ROOT  TO  THE  KNEW  ARRAY*  IF  L  «  0.  TRANSFER  THE  NSIMEC 
P  PARAMETERS  TO  THE  Q  ARRAY.  IF  L  «  1,  COMPOTE  A  NEW  SET  OF  S  NEC 

PARAMETERS  ANC  TRANSFER  THEM  1C  THE  Q  ARRAY.  SjNEC 

1  CO  2  1  =  1, N  SIMEC 

XNEWII»=XCLD(  i  )  NSINEC 

CO  2  J=i,N  NSINEC 

IF  (L.EQ.1I  GO  10  12  NSINEC 

Cl  I,JI*Pff,JI  NS  PEC 

GO  TO  2  NSINEC 

12  Cl  I,Jl=QCl  I,JMPtI,J)-QDtI,J)l*S/SF  NSINEC 

2  CONTINUE  NS IMEC 

CO  5  ITER*1,  ITMAX  NS IHEC 

— NS  IntG 

CALL  SUBROUTINE  JACOBI  TO  COMPUTE  THE  COEFFICIENTS  OF  THE  JACOBI  ANNS  IMEC 
MATRIX  ANC  STORE  THEM  IN  THE  A  ARRAY.  NSINEC 

- NSINEC 

CALL  JACGBI  (0, NMAX, XNEW,N,M* A )  NSINEC 

-  —  -•••--••••••---••••-•NS  Ir  EG 

CALL  SUBROUTINE  LSIMEQ  TO  COMPUTE  THE  SOLUTION  OF  A  SET  OF  NSINEC 

SIMULTANEOUS  LINEAR  ALGEBRAIC  EQUATIONS  USING  THE  GAUSS-JORCAN  NSINEC 
RECUCTION  SCHEME  WITH  THE  MAXIMUM  PIVOT  CRITERION.  THE  CONTENTS  NSINEC 
OF  THE  X  ARRAY  ARE  THE  CHANGES  WHICH  MUST  BE  NAOE  IN  THE  NSIMEC 

. ;2!!L«£!!: . ns!nec 

CALL  LSIMEQ  l  A, NMAX,  IR , JC,N,EPS1  ,X, IERR1  >  . . -NsImEC 

’  TEST'TO  CETERMINE  IF  THE  JACOBIAN  MATRIX  IS  SINGULAR  OR  NSINEC 

IL  LCONC  IT  ION  EC.  _  . . K|  IMEC 

”f~UERRULT~U  GC  TO  II  . . . . NSINEC 

. test~to~cItermInI~If  all  the  changes  in  thI  XNEW  ARRAY  ARE  NSINEC 

. !:!!!.!“!.!!S2; _ _ — . &f|S|g 

CO  3  1=1, N  NSINEC 

IF  I .NOT .ABSI XI I ) I .LT.EPS3)  GO  TO  A  NSIMEC 

3  CONTINUE  NSIMEC 

GO  TO  6  NSINEC 

- - - NSIMEC 

CONVERGENCE  HAS  NOT  BEEN  ACHIEVEO.  UPDATE  THE  XNEW  ARRAY  AND  NSINEC 
REPEAT  THE  ITERATION.  _  . . NSjMEC 

<,  CO  5  1  =  1, N  NSINEC 

5  xnewi  n=xNEwm«xm  M 5EIS 

MAXIMUM  NUMBER  OF  ITERATIONS  HAS  BEEN  EXCEECEO.  IF  THE  PARAMETER  NSINEC 
PERTURBATION  METHOD  IS  BEING  USEO,  HALVE  THE  STEP  SIZE  ANC  NSIMEQ 

CETERMINE  IF  IT  IS  LESS  THAN  EPS2.  . . N|  }MEC 

'"If"(lIIoIcT'gc"to"Io  ~~  NSIMEC 

S=S-CS  NSINEC 

CS  *0 . 5*CS  NSIMEC 

IF  (0S.LT.EPS2)  GC  TO  10  k  §  IMEC 

GO  TO  R  NSINEC 

- - - NSINEC 

CONVERGENCE  HAS  BEEN  ACHIEVEO.  IF  THE  PARAMETER  PERTURBATION  IS  NSINEC 

eEING  USEC,  TEST  TO  DETERMINE  IF  THE  STEPPING  PROCESS  IS  NSIMEC 

COMPLETE.  IF  IT  IS  OR  IF  THE  NEWTON  RAPHS0N  NETHOC  WAS  USEC,  NSIMEC 

RETURN  WITH  THE  SOLUTION  IN  THE  XNEW  ARRAY.  NSIMEC 

6  IF  IL.EQ.O)  GC  TO  9  NSIMEC 

IF  <  .NOT.S.LT.SF)  GO  TO  9  fcSIMEC 

TRANSFER  THE  NEW  ROOT  TO  THE  XCLO  ARRAY,  INCREASE  THE  STEP  SIZE  eYNSIMEC 
ITS  CURRENT  VALUE  AND  RETURN  TO  COMPUTE  NEW  C  PARAMETERS.  NSIMEQ 

--N5  I M EG 

CO  7  1=1, N  NSIMEC 

7  XOLCl  I  l-XNEWl  1 1  NSIMEC 

a  S-S40S  NSIMEC 

GO  TO  1  NSIMEC 

. . . NSINEQ 

RETURN.. .SUCCESSFUL  SOLUTION.  NS IMEC 

—  —  IrEG 

9  I6PP2=*1  NSIMEC 

RETURN  NSINEC 

r - - - NSINEC 
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C  ERROR  RETURN. ..IF  THE  NEWTON  RAPHSON  METHOO  HAS  iNTW| ESf RIfil IT= B  MsiSiS 

l  Si»X?'<8iiE8e?SoJ,SS,i^  "Slo,  IRe'IkI  llh  I?  HfalEIMBtg  ffl 

-  nsimec 

10  IERR2»-1  NS  I M EC 

11  RETURN  KsiSie 

C  *mm*«««m»m«**m«*m«m**;*****J*******5***mm*m*mmm”  ,r..EC 

SUBROUTINE  LSIHEQ  (R.NNAX,IR,JC,N,EPS1»X.IERR1I  LSIREC 

\ . i^r.i^misf^iSEH^rSsfsi^i^s^l-^^srHESjiftSE005  tug} 

C  SCHEME  WITH  THE  MAXIMUM_Pl VCT_CR I TERI ON. _ LSIHEQ 

C  DIMENSION  aTnmIxi 6 K IrInT » JcInT »xTn^ _ ISIMEC 

C  BEG In”thE  IElImINATION  PROCEOUrI.  _ LsJmEC 

e““Rwir: -  ttisa 

. _ LSIMEC 

c . BEgIn  THe” SEARCH  FOR  THE  MAXIMUM  PIVOT  ELEMENT.  _ LSIMEC 

c - -  lsimec 

BICA»0.0  LSIMEC 

CO  3  I*1#N  LSIMEC 

CO  3  J=1»N  lsimec 

IF  (KMl.Ic.O)  GO  TC  2 _ LSIMEC 

C  CHECK  ROw”aNC*COLUMN  PIVOT  SUBSCRIPTS  ALREADY  USED. _ LsImEC 

C . ------------T -  LSIMEC 

Jf  1  a  U(  in  go  TO  3  lsImec 

IF  (j'lEQ.klJjn  GO  TO  3  LSIMEC 

\  IF^( JnOT.ABSI  A( I* J 1 1 .GT.BIGAI  GC  TC  3  L S I M EG 

BIGA-ABSI A!  It Jll  LSIMEC 

1RIKI  *  LSIMEC 

JC  IK ;  *■  )  LSIMEC 

3  CONTINUE  LSIMEC 

C . CHicrT0'siriF’THE’MAxiMUM_pivCT_ELEMEN7_IS_GREATER_THAN_EPSl. _ LSjMEQ 

C  I F *7 TnOT T elG aTltII PSI*»  GO  TC  A  LsImEQ 

IERR1=-1  LSIMEC 

^RETURN . LSIMEQ 

C  NORMAL liI”THE~P I VCT  ELEMENT.  . . LsImEC 

C - - - — - — — - — —  LSIMEC 

4  IRK=  IRIK  I  LS  MEG 

JCK=JC  K»  LSIMEC 

eiGA»A(IRK,JCKI  LSIMEC 

5  S(IRKij|S!(i«K,JI/M6» _ tliSIS 

c  II iUiSJti'tHriiS  uu  IleSIkTs  In  thi  jciri  lsIheo 

g. ...... — -  LSIMEC 

IF  (i!eqIiRK{  GO  TO  7  LSIMEC 

AJCK’MI.JOO  LSIHEC 

6  AIUJl’illJh-AJCMAlIRK.JI  tllSlS 

7  CONTINUE _ _ _ LSIMEQ 

C  R  EC  RC)lR_THi_  SOLUTION  AND  TRANSFER  IT  TO  THE  X_ARR*Y. _ _ LSIMEC 

C - - - - - — — — — ~  ”  ”  “  LSIMEQ 

COS  I=1.N  LSIMEC 

IR I- IRI I )  LSIMEC 

8  XUCI  I* A<  IRI.NP1I _ LSIM^S 

c . Iuccessful’returnT  . lI  in  18 

C - - -  LSIMEC 

lERRlf *1  LSIMEC 

RETURN  LSIMEC 

SUBROUTINE  GVCINTIX.FXfXT.FTf NP,NOI  CVCINT 

CIMENSION  XTINPI.FTINPlf Ti 161  CVCINT 

N«NC 
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31  Nl«(N-l)/2 
N2=N/2 
N3  =NP-N2  + 1 
IF ( NP-NI 3C*  41*41 

41  N4=N1*2 

IF(Xfm-XT<  2)  122,80,60 
22  CONTINUE 

IFJX-2.*XTClHXT(2M20f20,2: 

21  IF(X-2.*XT(NP)4XT(NP-in441,441,20 
<•41  IFINP.LT.  10 ICO  TO  42 
N5*NP-N 
<*43  N5-N5/2 
N6-N44N5 

IF  (XKN6I.LT  ,XIN4*N6 
IF(N5.GT .1 ) GC  TO  443 

42  IF(X-XTIN4il45f43,43 

43  IF(N4-N3>44, 45*44 

44  N4*N!4*1 
GOTO  42 

45  N4rNW 
N5  =  N4-N1 
LJ46  1  =  1. N 

T (  1 1  =FT(N5 ) 

*6  N5=N54l 
L=  ( N*  1 1/  2 
TR*TIL  I 
N6=N4 
N7  =N44 1 
JU  =  1 
N2=N-1 
UN= 1 .0 
CO  12 J  =  1, N2 
N5*N4-N1 
N3=N-J 
COS  I  =  1,N3 

N8*N54J 

T(I)  =  (T(  Kl)-T(in/(XT(N8I-XT(N5»I 
9  ,\5=N54l 

GOTO!  lOtlll.JU 
1J  UN=UNM  X-XTi N6 I ) 

JU-2 
N6=N6-i 
GOTO  \2 

11  UN=UN*IX-XT(N7)) 

JU  =  1 

N7=N7-H 
L  =  L- 1 

12  TR=TR*UN4T(L I 
FX  =  TR 
RETURN 

20  MR  IT  El  6,  50  )  X,XTm,XT(NP» 

STOP 

50  FORMAT ( 23H  ARC.  NOT  IN  TABLE  X*  .E14.7.9H  XT(1)=  * 

1  E14. 7,101'  XT(NP I  *  ,E14.7*2X,6HDVDINT) 

20  MR  IT  E ( 6* 5 1 1  NP,ND 

51  FORMATI22H  TABLE  TOO  SHALL  NP»  ,I5»6H  ND*  ,I5,2X,6HCVCINT) 
STOP 

60  IF(X-2.*XTm4XT(2»»61, 20.20 

61  IF (X-2.*XT(NP  HXTINP-l 1  )20*721*721 
721  IFINP.LT. 10 ) CC  TO  72 

N5=NP-N 
723  N5=N5/2 
N6=N44N5 

IF(XT(N6I.GT.X)N4*N6 
IFIN5.GT.DGC  TC  723 

72  IFIX-XKN4)  173,73,45 

73  IF(N4-N3)74,45,74 

74  N4=N44l 
GOT 3  72 

80  MR  IT  E( 6*  52 1  XT(1) 

STOP 

52  FORMAT! 23F  CONSTANT  TABLE  XT(1)=  *E14.7,2X,6HDV0INT) 

ENC 

*4**4********4**********»* ******* «**«** *************  *******  *444*1 

SUBROUTINE  SCRTXY(X,Y,N> 

C I  HENS  ION  X ( N  ) , Y ( N ) 

R«N 

1  N*M/2 

IF  I H  .  EQ.G (RETURN 

K=N-M4  1 

J=1 

2  !  =  J 


>721,721 


CVCINT  4 

CVCINT  5 

CVCINT  fc 

CVCINT  7 

CVCINT  8 

CVCINT  S 

CVCINT1C 

CVCINT11 

CVCINT12 

CVCINT12 

CVCINT14 

CVCINT15 

CVDINT16 

CVCINT17 

CVCINT1B 

CVO  INTIS 

CVCINT20 

CVCINT21 

CVCINT22 

CVCINT23 

CVCINT24 

CVCINT25 

CVCINT26 

CVCINT27 

CVCINT28 

CVCINT2S 

CVCINT30 

CVCINT31 

CVCINT32 

CVCINT33 

CVCINT34 

CVCINT35 

CVCINT36 

CVCINT37 

CVCINT38 

CVCINT3S 

CVCINT40 

CVCINT41 

CVCINT4* 

CVCINT43 

CVCINT44 

CVCINT45 

CVCINT46 

CVCINT47 

CVCINT48 

CVCINT4S 

CVCINT50 

CVCINT51 

CVCINT52 

CVCINT52 

CVCINT54 

CVCINT55 

CVCINT56 

CVCINT57 

CVCINT58 

CVCINT59 

CVCINT60 

CVCINT61 

CVCINT62 

CVCINT63 

CVCINT64 

CVCINT65 

CVCINT66 

CVCINT67 

CVCINT68 

CVCINT69 

CVCINT70 

CVCINT71 

CVCINT72 

CVCINT73 

CVCINT74 

CVCINT75 

CVCINT76 

*•4*  1 

SORTXY  2 
SORTXY  3 
SORTXY  4 
SORTXY  5 
SORTXY  6 
SORTXY  7 
SORTXY  8 
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•  *VWWW,*W 


L-I+M 
Iff 


(XI  II.GT.XtUIGOTO  5 

4  J»J*1 

5  «{«’2,1,‘ 

5111:?"’ 

T-Y(LI 

Y(L)*Y(II 

INTERCHANGE  THE  I-TH  AND  L-TH  ELEMENTS  OF  ADDITIONAL  VECTORS. 
6  I-I-M 

IF( 1)4,4, 2 
ENC 


FUNCTION  TDINTiX.Y.XT.NX.YT.NV.FT.NXNAX.NFX.NPY) 
Cl  MENS  ION  XT(NX) ,YT(NY I ,rT(NXMAX»NY) ,F (161 «G(16) 


NX 1  =  NX 
NY  1  =  NY 
NPX1=NPX 
NPY1=NPY 
X1*X 
Y1=Y 

IFINPXi.GT .NX1.0R.NPY1.GT.NY1 IGCTO  14 
A*  2 .*XT( 1 1— XT 121 
E*  2 •*XT(NX1 1-XT(NXl-l) 

C*2.*YT(  1I-YK2J 
C*2.*YT(NY1)-YT(NY1-1) 

IF (XT( 1 ) «GT •XT(2) IGQTO  7 
IF(X1.LT.A.OR.X1.GT.B)GOTC  20 
CO  1  1=1, NX1 
IF (X 1 .LE  .XT( 1 1 IGOTQ  2 

1  CONTINUE 

2  I*  I-NPX1/ 2 

IFd.LT. 1)1*1  .  ,  , 

IFU.GT.NX1-NPX1*1H*NX1-NPX1*1 
IF(YT(1».GT.YT(2»)G0T0  9 
IFIY1.LT .C.0R.V1.GT.0IG0TC  21 
CO  2  J=1,NY1 
IF(Y1.LE.YT(JI)GOTO  4 

3  CONTINUE 

4  J* J-NPY 1/2 
IF ( J.LT • 1  IJ  =  1 

IF(J.GT.NY1-NPY1*1IJ*NY1-NPY1-»1 
IF(NPX1.LT .NPY1 IGCTO  11 
CO  6  K= 1 , NPY 1 
J1=J+K-1 
CO  5  L=1,NPX1 
11  =  I  ♦L- 1 

5  G( L)  =  FT(  IltJll 

6  CALL  " 

CALL 
RETURN 

7  IF(X1.GT.A.0R.X1.LT.BIG0TC  20 
CO  8  1=1, NX1 
IF(X1.GE.XT(  DIGOTG  2 

8  CONTINUE 
GOTO  2 

9  IFUl.GT.C.OR.Yl.LT.OIGOTC  21 
CO  10  J= 1 ,NY 1 
IF ( Y 1 .GE  . Y T ( J I (GOTO  4 

1G  CONTINUE 
GOTO  4 

11  CO  13  L= 1»NPX1 
I1=I*L-1 
CO  12  K= 1,NPY1 
J1=J-*K-1 

12  G(K)  =  FT(  U,J1I 

13  CALL 
CALL  _ 

RETURN 

20  MRITEI6, 221X1, XT(lt,XT(NXl) 

21  WR<{tE(6,22IY1,YT(1I,VT(NY1I 
STOP 

22  FOPMAT(30F  ERRCR  TDINT  X  Ct I  CF  TABLE  X«,E15.7,7H  XT (1I»,E15.7 ,8F 

231F0RMAT(l0H5ERR0R  TDINT  Y  CUI  CF  TABLE  Y=,E15.7,7H  YT (1) «,C15.7,8H 
1YT (NY )  =  , E 15 .7  ) 

14  KRITE(6,15INX1,NPX1,NY1,NPY1 

15  FORMAT! 32P  ERRCR  TDINT  TABLE  TCG  SMALL  NX=,I5»5H  NPX*,I5,4H  NY* 
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CCINT^(Y1,F(L I,YT, J,G, l,NPtl,NPYl,l,l) 
CC INTP(X1, TOINT,XT, I,r  ,l,NPXl,NPXl,l,ll 


SORTXY  9 

Iqrtxyiq 

S0RTXY11 
SQRTXYlf 
SORTXYl? 
S0RTXY14 
S9RTXY15 
S0RTXY16 
SORTXY 17 
S0RTXY1E 
S0RTXY19 
S0RTXY20 
S0RTXY21 
SORTXY22 


TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
TCINT 
, I5TCINT 


1 

2 

■a 

4 

c 

6 

7 

8 
9 

1C 

11 

12 

12 

14 

15 

16 

17 

18 
19 
2G 
21 
22 

23 

24 

25 

26 

27 

28 
29 
3C 

31 

32 
32 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


1Inc 


NPY*, 151 


TC1NT 

TCINT 


6e 

6S 


SUBROUTINE  OCINTPUX.FX.XT.NX.FT.NF.NPS.NC 
C IMENS ION  XT!NPS)»FTtNPS),Ml6l»S<lo) 


,1X1  *  IF  I ) 


41 


cl 


^43 


444 

42 

43 

44 


45 


x=xx 

NSX=NX 
NSF=NF 
NP=NPS 
N*NC 
IX- IX  I 
IV  - 1 F I 
Nl=(N-l)/2 
N2«N/2 
N3«NP-N2*1 
IFINP-NI30, 41,41 
N4=N 1*2 

N44=(N4-1)*IX*NSX 
N20=NSX* IX 
N22=C  NP-1 )*  IX+NSX 

*IF(XT(NSX*-XTIN20n22«80,|C 
IF <X-2.*XT(N$X)*XT(N20  20,^0,21 

IF (X-2.*XT!N22)4XTIN21)) 441, 441,20 
IF  ( NP  .LT • 10  )  GOTO  42 
N5  *NP-N 
N5=N5/2 
N6=N4+N5 

N66-I N6-1 1*  IX4NSX 
IF(XT(N66).GE.X)GCT0  444 
N4=N6 
N44=N66 

IFIN5.GT .2 IGCTO  443 
IFIX-XT1N44) 145,43,43 
IF!  N4-N 3)44, 45,44 
N4=N44l 
N4A=N44-»IX 
GOTO  42 
N4=N4~1 
N5  =N4-N1 

N20=(N5-1)*IY4NSF 
N21=(N5-1)*IX4NSX 
CO  46  I- l.N 
T (  I )  =  FT( N20 ) 

S(  I)=XTIN21) 

N2C=N20* IY 
N21=N21*IX 
1=  !  N+ 1 )/ 2 
TR=T1L ) 

N4=N  1*1 
N6=N4 
N7=N4«1 
JU  =  i 
N2=N- 1 

UN=1.G  „ 

IF ( N • EO , 1 ) GOTO  13 
CO  12  J=1.N2 
N5=N4-N1 
N3=N- J 
CO  9  1=1, N3 

T(  n  =  (T(  1*1  l-T<  I))/tSIN8)-SIN5)) 

N5=N5*1 

GOTO! 10, 11), JU 
UN=UN*IX-S!N6 ) ) 

JU  =  2 
N6=N6-1 
GOTO  12 

UN=UN*(X-S(N7I) 

JU  =  1 
N7=N7«1 
L-L-l  ,  . 

TR«TR*UN*T(L) 

FX  =  TR 
RETURN 

WRITE(6,50IX,XTINSX),XT(N22) 

50  F0RMATI23F  ARG.  NOT  IN  TABLE  X*  |Ei5?7*9M 
1E14.7.10F  XT (NPS  )  =  ,E14, 7, 2X, 6H00INTP ) 

51  FORMAT! *Ta6lE  TOO  SMALL  NPS=  ,15, 6H  ND=  , 1 5 ,2X ,6HDD INT F ) 

STOP  82 


46 


10 


11 


12 

13 

20 


XT  ( 1 )  = 


t 

i 

4 

c 

6 

7 

8 
<3 


***♦ 

CCINTP 
CC1NTP 
CCINTP 
CCINTP 
CCINTP 
CCINTP 
CCINTP 
CCINTP  . 
CCINTP10 
CCINTP11 
CC  INTP12 
CCINTP13 
CCINTP14 
CCINTP15 
CCINTP16 
CCINTP17 
CCINTP18 
CCINTPIS 
rCINTP2C 
CCINTP21 
CCINTP22 
CCINTP22 
CCINTP24 
CC INTP25 
CCINTP26 
CCINTP27 
CL INTP28 
CCINTP2S 
CC INTP3C 
CCINTP31 
CCINTP32 
CCINTP33 
CCINTP34 
CCINTP35 
CCINTP36 
CC INTP37 
LL INTP38 
CCINTP39 
CCINTP40 
CCINTP41 
CCINTP42 
CCINTP43 
CCINTP44 
CC1NTP45 
CCINTP46 
CCINTP47 
CCINTP4P 
CCINTP4S 
CCINTP50 
CCINTP51 
CCINTP52 
CCINTP53 
CCINTP54 
CCINTP55 
CCINTP56 
CCINTP57 
CCINTP5E 
CCINTP5S 
CCINTP60 
CCINTP61 
CCINTP62 
CCINTP62 
CCINTP64 
CCINTP65 
CCINTP66 
CCINTP67 
CCINTP6E 
CCINTP6S 
CCINTP7C 
CCINTP71 
CCINTP72 
CCINTP73 
CCINTP74 
CCINTP75 
CCINTP76 
CCINTP77 
CCINTP76 
CCINTP7S 


IF (X-2.*XT(NSX)4XT(N2C  61,20, 2C 
IF(X-2. *XT (N 22I ♦XT(N21II 20. 721.721 
IF  CNP.LT.IOIGOTC  72 


N5-NP-N 

723  N5»N5/2 
N6=N4*N5 

N66=(N6-1I*IX4NSX 
IF(XT(N66I.LE.X)GCTC  724 
N4  =  N6 
N44=N66 

724  IF IN5.GT .2 1G0TC  723 

72  IF ( X-XT ( N44 H73, 73.45 

73  IFIN4-N3I74.45.74 

74  N4=N4*1 
N44=N44*  IX 
GOTO  72 

cO  *R  ITE(6,52IXT(NSX J 

•=2  FORMAT ( 2  3F  CONSTANT  TABLE  XT<  1 1  =  ,E  14.7 ,2X .6H00I MP) 

^SiisS^Msa^SiiiiiaSiSiSi™ 
,inSWKi®ami:sm . . 

ITTC=0 

INCC-0 

I F (  ITYP.GT.1IGCTO  7C 
ERM=5 .*ER 
ER 1= . 1*ERM 
CPST-CNXT 
1=1 

IFf IrIIeG.O. »ER2=.00C001*PS 

10  GOTO  (  20i90|[llC»  13C,  150.7C.25C1.I 
20  IF  t ITYP.Ic.l IREtORN 
25  CALL  TERMY.YP.TC.PVI 
GOTO  (30,2701.1 

31  SStt  KIM-ici,, 

IF ( TCI IiIStIc.  IGCTC  41 
FB ( I  Is 1 • 

F  A (  11  =  0. 

GOTO  SO 
41  F  A (  11  =  1, 

FB (  11=0. 

50  CONTINUE 

PV2=AINT(PV/PS*SIGN(ER2,ONXT)  I 

IF ( ABSCPV2-PV I .LT.ER2 JPV2*PV2-PS 
K*  0 

70  F6=CNXT/6. 

F3=F6-»H6 
F8=.125*CNXT 
F2=.5*CNXT 
CO  80  1=1. N 
UII=Y(II 

ci (  i»=yp(  1  j 

eo  v«  n=c<  n«N3*Qi(  i » 

1  =  2 

goto  10 

*C0  Y? I i  =  0( I >  *H6*  €  QIC  I ) ♦ YP 1 1  1 1 
1=3 

no  coTi20  i=i.n 

12c  Y(l!=Q?n*HeiiQ2m*Qim» 

1*4 

GOTO  10 

130  f?4i?iUo|j|! 

140  Y(  ii=C(U«H2*(Tl«Q2CI)» 

1=5 

GOTO  10 
150  ERHAX=0. 

CO  19C  1=1. N 


CCINTP8C 
CCINTPB1 
CCINTP82 
CCINTP63 
CCINTP84 
CCINTP85 
CCINTP86 
CCINTP87 
CCINTP88 
CCINTP88 
CCINTP9C 
CCINTP91 
CCINTP92 
CC1NTP92 
CCINTP94 
CCINTP95 
CC  INTP9c 
C  C  I N  (  P  9  7 
CCINTPTE 
CC  INTP99 


KUTHF  1 
KUTHE  2 
kutme  3 

KUTHE  4 
KUTHF  5 
KUTHE  t 
KUTHE  7 
KUTHE  e 
KLi  THE  9 
KUTHE  10 
KUTHE  11 
KUTHE  12 
KUTHE  12 
KUTHE  14 
KUTPF  15 
KUTHE  U 
kuthe  17 
kuthe  if 

KUTHE  18 
KUTHF  2C 
KUTHE  21 
KUTHE  22 
KUTPF  23 
KUTHE  24 
KUTHE  25 
KUTHE  26 
KUTME  27 
KUTME  26 
KUTME  28 

KUTME  31 
KUTME  32 
KUTMFA32 
KUTHE  33 
KUTHE  34 
KUTHE  35 
KUTHE  36 
KUTHE  37 
KUTHE  38 
KUTME  39 
KUTHE  40 
KUTHE  41 
KUTHE  42 
KUTHE  42 
KUTHE  44 
KUTHE  45 
KUTHE  46 
KUTHE  47 
KUTHE  48 
KUTHE  48 
KUTHE  50 
KUTHE  51 
KUTHE  52 
KUTHE  52 
KUTHE  54 
KUTHE  55 
KUTHE  56 
KUTHE  57 
KUTHE  56 
KUTHE  58 
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KUTRE  fcC 


1 1  =ci  n«Q2in«YP(m 

T ,  -  r  A  B  f.  (  1 1  -  Y  I  I  I  ) 

'  <  I  I  =  T  I 

I  f  (  RM.EC.G.  IGCTO  19C 
CO  TO  I  170, 16C,  lFCl.MF 

0  IF(ABS(T1I.GE.1.»T2=I2/*BSI III 
.70  If(T<..GT.ERMAX)ERHAX=T2 

GOTO  Is 0 

.60  IF1T1.EQ.C.  IGCTO  170 
T i  =  T 2/CAESI T 1 I 
l-H 

I2=WC  J-1MT2 
Ti*W(  J)M1 
IF( T2.CT.T1) I2=T1 
GOTO  170 
nO  rONT  INllE 
tPST -CNXT 

I F (  EKM.fc.C  •  »GCTC  29C 

II  (  ITTC.EC.l )G0TC  2m0 
FT (cRRAX.EQ.C.ICCT'J  21C 

I  1  =FR1/[R"AX 

ll-ITl.GI  . i . .ANG. INDD.GT.CIGCTC  192 

I I  ( «RS( T 1-1. 1. 1 T.. 5 IGCTO  195 
i  V<T-CPGl*T  1  **  •  2 

•  'H  ;  2 1  i 
.  IN  C  l‘  =  1M0I>  1 

'  T'  2dj 

‘  .  :  \  X  i  -  C  P  S  T  ♦  (  1 .  ♦  .  2*  I  fl-l.H 
•.  IF ( "HSIDNXT I .LT.CRAX  IGCTC  2^0 
If  (  INCC.GT.u IGCTG  192 
C.itdr.E'AX 

iriOPbT.LT.C . II  NX  T --DM AX 

if  (ekpax.lt. t«M igcto  29c 

GO  :2C  I  I  *  N 
.3,.  i  I  Ihun 
inrrd 
1  -  c 

GOTO  10 
'  9  ,  !  --  7 

ih  ityp.ci.iireturn 

TPV=PV 

CO  260  l=l.NTC 
.  6C  .11  I  )dC(  I) 

1  =  2 

COT'  25 

,  7*.  CO  2oO  I  -  1, NTC 

I  F  ( F  A  (  H.EQ.FEI  I)  IGCTO  37C 
IFITCI  M.GT.C.CC)  GC  TC  28C 
r  b  C  11  =  1. 

IFIABSNim  I.GE.CR2.0R.K.GT.  1 IGCIC  290 
F  A { II  =  0 . 

COT C  290 
-EG  F A (  l|=i. 

IF ( APSIQK  I  I I. GE. ER2. OR. K.GT.l IGCTC  290 
FBI  11=0. 

.•>U  IFIFA4  II.EQ.FBI  IllGOTO  37C 

3CC  CONTINUE 
K  =  K  ♦  1 

!F(l'eS(PV-PVl).lT.ER2)GClC  j5C 
IF ( A6SIPV-PV2 I.LT .ER  2 IGCTC  25C 
IFire. NF.;. IGCTO  290 
;G5  IF1FV.GT.PV1  I0CT0  330 
IF (PV.GT.PV2 IGCTG  70 
PV  2  PV2 
-.10  F  G  -  1  . 

i,n  =  r.f!XT 

j 2G  CNXT-CPSTMPV2-PV l/IPV-TPVI 
IT  TC  =  1 

GOTO  70 
2 20  PV  2 -PV 1 
con  310 

^90  IFIFC.GT.9. IGCTC  350 
FC=FC*1. 

COT''  320 
j  50  1  =  2 

GOTO  30 
it 0  PV1  =  PV4Pj 
FV  2  =PV-P 9 
I T  T  C  =  0 

IF( FC.EQ.C. I< C TC  2C5  , 

IC=0.  1 
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KUIRE  c< 
KUTRF  62 
KUTKE  69 
K  U  T  B  E  65 
KUTKE  66 
KUTKE  67 
KUTRE  66 

KUTRE  7C 
KUTBF  71 
KUTRE  7  l 
Kill  PE  7  2 
K  LITRE  79 
KUTRE  7 5 
KUTRF  It 
KUTRE  77 
KUTRE  78 
KUTRE  79 
KITRE  8C 
KLlt-E  8  1 
KUTRF  82 
p. I . T M  F  82 
KUTRE  89 
KUTRE  85 
KUTRE  8c 
KOTMt  87 
KUTRE  86 
KUTRF  89 
KUTRE  9C 
KUTRE  91 
KlJ  T 8 E  92 
KUTRE  92 
KUTRF  99 
KUTRF  95 
KUTRF  96 
KUTRE  97 
KUTRF  96 
KJTRE  98 
HUTRE10C 
HUT  PE  10  1 
KtJTHE  10  2 
HUTRE102 
KUIFE109 
KUTRE105 
KLTRE106 
HUT  8  E 107 
KUTRE106 
KtIFn09 
KUTME110 
KUTRE111 
KUTHE112 
HUT8F113 
KUTRE119 
KUTRF115 
KUTRE116 
KUTRE117 
KUTMA117 
KUTRE118 
KUTRE119 
KIJTME12G 
KUTPE121 
KUTRE122 
KLTME122 
KUTRE129 
KUTRE125 
KUTRE126 
KITRE127 
KUT8E128 
KUTRE129 
KUTME13C 
KUTRE131 
KUTRE132 
KUTBE133 
KUTME139 
KUTRE135 
KUTRE136 
KUTRE137 
hUT8E13e 
KUTRE139 
KUTRE19C 


I 


CNXT-SH 
GOTO  305 

370  IF(«eS(TCUi  I.LT.ER2IG0T0  360 
cnxt*cpswc<  n/coim-Tcun 

IT  T  C  *  1 
GOTO  70 
380  NT S=  I 
1*3 

GOTO  30 
390  RETURN 
ENC 

J .15  1.66 

20.0 


u  .40 
G  .25 
0.2 


1.76 

1.66 

1.66 


KUTME141 

KUTPE141 

KUTME143 

KUTPE144 

KUTNE145 

KUTPE146 

KUTME147 

KUTPEl^e 

KUTPE149 

KUTPE15G 

KUTME151 


1.43 

.159 

.349 

.302 

-.385 

0.04 

-1.4 

.159 

.349 

.302 

-.385 

0.04 

-1.4 

.159 

.349 

.302 

-.385 

0.04 

1.45 

.159 

.349 

.302 

-.385 

0.04 

85 


« 


LIST  OF  SYMBOLS 


a  radius  of  circular  cylinder  in  physical  plane 

2 

Cf.  skin  friction  coefficient,  (3u/3y)w/ (puQ  /2) 

d  diameter  of  rotating  cylinder  in  physical  plane 

h  coordinate  scaling  function  of  x,  n  =  y/h(x) 

i  unit  imaginary  value 

k  nondimens ionali zed  outer  velocity  at  the  separation  point  on 

the  upper  part  of  the  cylinder 

m  angle  used  in  z  plane,  see  Figure  2 

mg  angle  of  m  corresponding  to  separation 

p  local  pressure 

q  inviscid  velocity  along  cylinder  surface  in  f;  plane 

q  dimensioned  inviscid  velocity 

r,<*>  polar  coordinates  in  I  plane 

r ,o  polar  coordinates  in  1  plane 

u  nondimensionaiizsd  local  velocity  in  direction  tangential  to 

cylinder 

u  nondimensionalized  circulation  velocity  on  circle  in  £  plane 

c 

u  inviscid  velocity  along  cylinder  surface  in  z  plane 

uw  peripheral  velocity  of  cylinder  wall 

u  free-stream  velocity 

o 

u  dimensioned  form  of  u 

e  e 

0  dimensioned  boundary- layer  velocity  in  direction  tangential 

to  cylinder 

v  boundary-layer  velocity  in  direction  normal  to  cylinder 

dimensioned  local  velocity  in  direction  normal  to  cylinder 
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LIST  OF  SYMBOLS  (Continued) 
wf )  complex  inviscid  velocity  in  plane 

x  nondimensional ized  coordinate  along  surface  of  cylinder  -  X/n 

%  distance  along  cylinder  from  X  axis 

v  nondimens iona I i zed  coordinate  normal  to  surface  -  origin  at 

surface  of  cylinder  in  physical  plane 

y  dimensioned  value  of  y 

complex  coordinates  in  physical  plane  z  =  X  +  i  Y 

complex  coordinate  rotated  by  angle  -•  to  the  ^  plane, 

'  =  X  +  i  Y 

UU  flu  drag  coefficient  on  the  cylinder 

C,  the  lift  coefficient  on  the  cylinder 

n 

F  value  of  the  convergence  criterion 

c 

llj  asymptotic  breadth  of  wake 

Ml”)  conformal  mapping  function,  maps  from  '  plane  to  z  plane 

N;:j  function  relating  the  fluid  velocities  on  the  untransformed 

and  transformed  circles 

Q  nondimensional ized  strength  of  single  source  on  upper  part 

of  cylinder 

O  (  nondimensionali zed  strength  of  single  source  on  lower  part 

of  cy 1 i nder 

0 j  dimensioned  value  of 

U  ,  1  i  mens  i  oned  va  1  ue  o  f  Q 

d  radius  of  circle  in  plane,  R  =  esc  > 

S.  upper  separation  point 

X 

S  ,  lower  separation  point 
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LIST  OF  SYMBOLS  (Continued) 


magnitude  of  free-stream  velocity  in  5  plane 

complex  velocity  in  z  plane 

polar  angle  of  upper  separation  point  in  \  plane 

angle  of  upper  source  from  real  axis  in  5  plane 

angle  of  lower  source  from  real  axis  in  c  plane 

angle  that  \  system  needs  to  be  turned  to  coincide  with  ; 
system 

angle  that  z  system  needs  to  be  turned  to  coincide  with  z 
system  after  a  translation  of  z  coordinates  so  that  the 
origins  of  the  two  systems  coincide 

primary  untransformed  complex  plane 

untransformed  complex  plane  rotated  by  an  angle  of  -e  to 
the  c,  plane 

viscosity 

kinematic  viscosity 

real  coordinate  for  x,  plane 

real  coordinate  for  £  plane 

mass  per  unit  volume 
polar  angle  in  c  plane 
surface  shear  stress 
polar  angle  in  c  plane 
complex-potential  flow  field 
imaginary  coordinate  for  <;  plane 

imaginary  coordinate  for  plane 


circulation  strength 
flow  potential 
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I  1ST  OF  SYMBOLS  (Continued) 
angular  increment  from  separation  points 
SUBSCRIPTS 

e  boundary-layer  edge  conditions 

f  final  conditions,  also  stands  for  fixed  wall  conditions 

i  initial  conditions 

s  boundary- layer  separation  conditions 

u  wall  conditions 

1  upper  side  conditions 

~  lower  side  conditions 

SUPliRSCRlPTS 

some  dimensioned  quantities 

transformed  quantities  with  a  few  exce|tions 
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